4 Electromagnetic wave propagation in dielectrics 



4.1 Introduction 

It is easily demonstrated that the fields associated with an electromagnetic wave 
propagating through a uniform dielectric medium of dielectric constant e satisfy 



e ^2 



- I = 0, (4.1) 



c2 dt'^ 



and 

VA-E = -^. (4.2) 

dt ^ ' 

The plane wave solutions to these equations are well known: 

E = Eoe'^^-''-''^\ (4.3a) 

B = Boe''^^-''-'^^\ (4.3b) 

where Eq and Bq are constant vectors, with 



a;2 c2 



A;2 e 

and 



a; 

The phase velocity of the wave is given by 

UJ c 
k n ' 



(4.4) 



Bo = (4.5) 



(4.6) 



where 

n = (4.7) 

is called the refractive index of the medium. It is clear that an electromagnetic 
wave propagates with a phase velocity which is slower than the velocity of light 
in a conventional {i.e., e real and greater than unity) dielectric medium. 
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In some dielectric media e is complex. This leads, from Eq. (4.4), to a complex 
wave vector k. For a wave propagating in the x-direction we obtain 

E = Eq exp[i (Re(A;) x - cut)] exp[-Im(A;) x]. (4.8) 

Thus, a complex dielectric constant leads to the attenuation (or amplification) of 
the wave as it propagates through the medium in question. 

Up to now, we have tacitly assumed that e is the same for waves of all frequen- 
cies. In practice, this is not the case. In dielectric media e is, in general, complex, 
and varies (in some cases, strongly) with the wave frequency, lu. Thus, waves of 
different frequencies propagate through a dielectric medium with different phase 
velocities. This phenomenon is known as dispersion. Moreover, there may exist 
frequency bands in which the waves are attenuated {i.e., absorbed). All of this 
makes the problem of determining the behaviour of a wave packet as it propagates 
through a dielectric medium far from straightforward. Recall, that the solution 
to this problem for a wave packet traveling through a vacuum is fairly trivial. 
The packet propagates at the velocity c without changing its shape. What is the 
equivalent result for the case of a dielectric medium? This is an important ques- 
tion, since nearly all of our information regarding the universe is obtained from 
the study of electromagnetic waves emitted by distant objects. All of these waves 
have to propagate through dispersive media {e.g., the interstellar medium, the 
ionosphere, the atmosphere) before reaching us. It is, therefore, vitally important 
that we understand which aspects of these wave signals are predominantly de- 
termined by the wave sources, and which are strongly modified by the dispersive 
media through which they have propagated in order to reach us. 

The study of wave propagation through dispersive media was pioneered by 
two scientists, Arnold Sommerfeld and Leon Brillouin, during the first half of 
this century. In the following discussion, we shall stick as close as possible to 
Sommerfeld and Brillouin's original analysis. 

4.2 The form of the dielectric constant 

Let us investigate an electromagnetic wave propagating through a transparent, 
isotropic, non-conducting, medium. The electric displacement inside the medium 
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is given by 

D = eoE + P, (4.9) 

where P is the electric polarization. Since electrons are much lighter than ions 
(or atomic nuclei), we would expect the former to displace further than the lat- 
ter under the influence of an electric field. Thus, to a first approximation the 
polarization P is determined by the electron response to the wave. Suppose that 
the electrons displace a distance s from their rest positions in the presence of the 
wave. It follows that 

P = -Nes, (4.10) 
where N is the number density of electrons. 

Let us assume that the electrons are bound "quasi-elastically" to their rest 
positions, so that they seek to return to these positions when displaced from them 
by a field E. It follows that s satisfies the differential equation of the form 

ms-\-fs = -eE, (4.11) 

where m is the electron mass, —fs is the restoring force, and ' denotes a partial 
derivative with respect to time. The above equation can also be written 

s + gcjQS + ujo^s = E, (4.12) 

m 

where 

/ 



Wo 



(4.13) 

m 

is the characteristic oscillation frequency of the electrons. In almost all dielec- 
tric media this frequency lies in the far ultraviolet region of the electromagnetic 
spectrum. In Eq. (4.12) we have added a phenomenological damping term gcuo s, 
in order to take into account the fact that an electron excited by an impulsive 
electric field does not oscillate for ever. In general, however, electrons in dielectric 
media can be regarded as high-Q oscillators, which is another way of saying that 
the dimensionless damping constant g is typically much less than unity. Thus, 
an electron "rings" for a long time after being excited by an impulse. 
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Let us assume that the electrons oscillate in sympathy with the wave at the 
wave frequency uj. It follows from Eq. (4.12) that 

« = -^M^ . (4.14) 

Note that we have neglected the response of the electrons to the magnetic compo- 
nent of the the wave. It is easily demonstrated that this is a good approximation 
provided that the electrons do not oscillate with relativistic velocities [i.e., pro- 
vided that the amplitude of the wave is sufficiently small). Thus, Eq. (4.10) 
yields 

P=^llfM^. (4.15) 

ujq — UJ''' —\gujujQ 



Since, by definition, 
it follows that 



D = tQtE = tQE ^ P, (4.16) 

uJq^ —uj^—\gujujQ 

Thus, the index of refraction is frequency dependent. Since ujq typically lies 
in the ultraviolet region of the spectrum (and since g <C 1), it is clear that 
the denominator ujq' — — igujujQ ~ ujq — uS^ is positive in the entire visible 
spectrum, and is larger at the red end than at the blue end. This implies that 
blue light is refracted more than red light. This is normal dispersion. Incidentally, 
an expression, like the above, which specifies the dispersion of waves propagating 
through some dielectric medium is usually called a dispersion relation. 

Let us now suppose that there are TV molecules per unit volume with Z elec- 
trons per molecule, and that instead of a single oscillation frequency for all elec- 
trons, there are fi electrons per molecule with oscillation frequency uji and damp- 
ing constant gi. It is easily demonstrated that 

n\u) = 1 + ^ V , (4.18) 

Com ^ - u'' -igiUUi 

where the oscillator strengths fi satisfy the sum rule, 

E/^ = ^- (4-19) 
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A more exact quantum mechanical treatment of the response of an atom, or 
molecule, to a low amplitude electromagnetic wave also leads to a dispersion 
relation of the form (4.18), except that the quantities /i, uji, and gi can, in 
principle, be calculated from first principles. In practice, this is too difficult 
except for the very simplest cases. 

Since the damping constants are generally small compared to unity, it 
follows from Eq. (4.18) that n(a;) is a predominately real quantity at most wave 
frequencies. The factor [co-'^ — co'^)~^ is positive for uj < coi and negative for 
uj > uJi- Thus, at low frequencies, below the smallest a;^, all of the terms in the 
sum in (4.18) are positive, and n{uS) is consequently greater than unity. As uo is 
raised so that it passes successive uJi values, more and more negative terms occur 
in the sum, until eventually the whole sum is negative and n{uS) is less than unity. 
Thus, at very high frequencies electromagnetic waves propagate through dielectric 
media with phase velocities which exceed the velocity of light in a vacuum. For 
UJ UJi, Eq. (4.18) predicts a rather violent variation of the refractive index with 
frequency. Let us examine this phenomenon more closely. 



4.3 Anomalous dispersion and resonant absorption 

When UJ is approximately equal to uJi the dispersion relation (4.18) reduces to 



■2 , Ne'^fi/eom 
n -\- — ^ 2 ^ (4-20) 



,2 „ 2 

uj.^ - uj^ - igiUJUJi 



where is the average contribution in the vicinity of u = Ui of all other res- 
onances (also included in is the contribution 1 of the vacuum displacement 
current, which was previously written down separately). The refractive index is 
clearly complex. For a wave propagating in the x-direction 

E = Eq exp[i {uj/c){Re{n) x - ct)] exp[-(a;/c) Im(n) x]. (4-21) 

Thus, the phase velocity of the wave is determined by the real part of the refractive 
index via 

^ = (4-22) 
Re(n) ^ ^ 
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Note that a positive imaginary component of the refractive index leads to the 
attenuation of the wave as it propagates. 



Let 



^2 ^ fi^ ^ (4.23a) 



Com ujj^ ^ ' 
a;2 - a;, 2 



X = (4.23b) 



UJ. 



a" 

2Re(n)Im(n) 
^2 



(4.23d) 



where a, x, y, z are all dimensionless quantities. It follows from Eq. (4.20) that 



'^^^ (4.24b) 



X2 + £/.2(l +x) 



Let us adopt the physical ordering <C 1. The extrema of the function y occur 
at a; = ztgi. It is easily demonstrated that 



n-2 1 

2/min = y(x = Qi) = -Y - (4.24c) 

n-2 1 

2/max = y{x = -9i) = ^ + 73— • (4.24d) 

a Zgi 

The maximum value of the function z occurs at a; = 0. In fact, 

1 

^max = — . (4.25) 
9i 

Note that 

z{x = ±gi) = ^_. (4.26) 
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-Si g. x-> 
Figure 5: Sketch of the variation of the functions y and z with x 



Figure 5 shows a sketch of the variation of the functions y and z with x. 
These curves are also indicative of the variation of Re(n) and Im(n), respectively, 
with frequency uj in the vicinity of the resonant frequency uji. Recall that nor- 
mal dispersion is associated with an increase in Re(n) with increasing u. The 
reverse situation is termed anomalous dispersion. It is clear from the figure that 
normal dispersion occurs everywhere except in the immediate neighbourhood of 
the resonant frequency uJi. It is also clear that the imaginary part of the refrac- 
tive index is only appreciable in those regions of the electromagnetic spectrum 
where anomalous dispersion takes place. A positive imaginary component of the 
refractive index implies that the wave is absorbed as it propagates through the 
medium, so the regions of the spectrum where Im(n) is appreciable are called 
regions of resonant absorption. Anomalous dispersion and resonant absorption 
take place in the vicinity of the ith resonance when \uj — uji\ < 0{gi). Since the 
damping constants Qi are, in practice, very small compared to unity, the regions 
of the spectrum in which resonant absorption takes place are strongly localized 
in the vicinity of the various resonant frequencies. 
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The dispersion relation (4.18) only takes electron resonances into account. 
Of course, there are also resonances associated with displacements of the ions 
(or atomic nuclei). The off-resonance contributions to the right-hand side of 
Eq. (4.18) from the ions are smaller than those from the electrons by a factor of 
order m/M (where M is a typical ion mass). Nevertheless, the ion contributions 
are important because they give rise to anomalous dispersion and resonant ab- 
sorption close to the ion resonant frequencies. The ion resonances associated with 
the stretching and bending of molecular bonds typically lie in the infrared region 
of the electromagnetic spectrum. Those associated with molecular rotation (these 
resonances only affect the dispersion relation if the molecule is polar) occur in 
the microwave region of the spectrum. Thus, both air and water exhibit strong 
resonant absorption of electromagnetic waves in both the ultraviolet and infrared 
regions of the spectrum. In the first case this is due to electron resonances, and 
in the second to ion resonances. The visible region of the spectrum exists as a 
narrow window lying between these two regions in which there is comparatively 
little attenuation of electromagnetic waves. 



4.4 Wave propagation through a conducting medium 

In the limit a; — ?► 0, there is a significant difference in the response of a dielectric 
medium, depending on whether the lowest resonant frequency is zero or non-zero. 
For insulators the lowest resonant frequency is different from zero. In this case, 
the low frequency refractive index is predominately real, and is also greater than 
unity. Suppose, however, that some fraction /q of the electrons are "free," in the 
sense of having ujq = 0. In this situation, the low frequency dielectric constant 
takes the form 

e{u:)^n\u,)=no' + i— , (4.27) 

Com a; (70 — 1 uj) 

where no is the contribution to the refractive index from all of the other reso- 
nances, and 7o = lim^^^Q gQCjQ. Note that for a conducting medium the contri- 
bution to the refractive index from the free electrons is singular at a; = 0. This 
singular behaviour can be explained as follows. Consider the Ampere-Maxwell 
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equation 

VAB = fio(jt + ^]. (4.28) 




Let us assume that the medium in question obeys Ohm's law, jt = aE, and has 
a "normal" dielectric constant Uq^. Here, a is the conductivity. Assuming an 
exp(— iwt) time dependence of all field quantities the above equation yields 

YAR = -ieo uj (rio^ + i — ) ^. (4.29) 

Suppose, however, that we do not explicitly use Ohm's law but, instead, attribute 
all of the properties of the medium to the dielectric constant. In this case, the 
effective dielectric constant of the medium is equivalent to the term in round 
brackets on the right-hand side of the above equation. Thus, 

e(u;) = n2(u;) = + i— . (4.30) 

Co CO 

A comparison of this term with Eq. (4.27) yields the following expression for the 
conductivity 

cr=— (4.31) 

m(7o - lo;) 

Thus, at low frequencies conductors possess predominately real conductivities 
(i.e., the current remains in phase with the electric field). However, at higher 
frequencies the conductivity becomes complex. At these sorts of frequencies there 
is little meaningful distinction between a conductor and an insulator, since the 
"conductivity" contribution to e(a;) appears as a resonant amplitude just like the 
other contributions. For a good conductor, such as Copper, the conductivity 
remains predominately real for all frequencies up to and including those in the 
microwave region of the electromagnetic spectrum. 

The conventional way in which to represent the complex refractive index of a 
conducting medium (in the low frequency limit) is to write it in terms of a real 
"normal" dielectric constant, e = n^^, and a real conductivity, a. Thus, from 
Eq. (4.30) 

n'^{uj) = e + i— . (4.32) 

Co CO 
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For a poor conductor (cr/e cq w <C 1) we find 

A; = n- ~ ^/^- + i— ^ . (4.33) 

c c 2^JeeQC 

In this limit Re(/c) Im(/c), and the attenuation of the wave, which is governed by 
lm.{k) [see Eq. (4.8)], is independent of the frequency. Thus, for a poor conductor 
the wave is basically the same as a wave propagating through a conventional 
dielectric with dielectric constant e, except that the wave attenuates gradually 
over a distance of very many wavelengths. For a good conductor {a/eeocu ^ 1) 

k ~ e^'^/^^A^^^. (4.34) 

It follows from Eq. (4.5) that 



^ = ^ = (4.35) 

Thus, the phase of the magnetic field lags that of the electric field by 45°. More- 
over, the magnitude of cBq is much larger than that of Eq (since a/tQUJ ^ e > 1). 
It follows that the field energy is almost entirely magnetic in nature. It is 
clear that an electromagnetic wave propagating through a good conductor has 
markedly different properties to a wave propagating through a conventional di- 
electric. For a wave propagating in the a^-direction, the amplitudes of the electric 
and magnetic fields attenuate like exp{—x/d), where 



d=J . (4.36) 

This quantity is known as the skin depth. It is clear that an electromagnetic wave 
incident on a conducting medium will not penetrate more than a few skin depths 
into that medium. 



4.5 The high frequency limit 

Consider the behaviour of the dispersion relation (4.18) in the high frequency 
limit UJ ^ uji (for all i). In this limit, the relation simplifies considerably to give 

n^H = 1 - ^, (4.37) 
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where the quantity 




NZe'^ 



Com 



(4.38) 



is called the plasma frequency. The wave-number in the high frequency limit is 
given by 



This expression is only valid in dielectrics when uj ujp. Thus, the refractive 
index is real and slightly less than unity, giving waves which propagate without 
attenuation with a phase velocity slightly larger than the velocity of light in vac- 
uum. However, in certain ionized media (in particular, in tenuous plasmas such 
as occur in the ionosphere) the electrons are free and the damping is negligible. 
In this case, Eqs. (4.37) and (4.39) are valid even when uj <ujp. It is clear that a 
wave can only propagate through a tenuous plasma if its frequency exceeds the 
plasma frequency (in which case it has a real wave-number). If wave frequency 
is less than the plasma frequency then the wave-number is purely imaginary, ac- 
cording to Eq. (4.39), and the wave is therefore attenuated. This accounts for the 
fact that long-wave and medium-wave radio signals can be received even when 
the transmitter lies over the horizon. The frequency of these waves is less than 
the plasma frequency of the ionosphere, which reflects them, so they are trapped 
between the ionosphere and the surface of the Earth (which is also a good re- 
flector of radio waves), and can, in certain cases, travel many times around the 
Earth before being attenuated. Unfortunately, this scheme does not work very 
well for medium- wave signals at night. The problem is that the plasma frequency 
of the ionosphere is proportional to the square root of the number density of free 
ionospheric electrons. These free electrons are generated through the ionization 
of neutral molecules by ultraviolet radiation from the Sun. Of course, there is no 
radiation from the Sun at night so the density of free electrons starts to drop as 
the electrons gradually recombine with ions in the ionosphere. Eventually, the 
plasma frequency of the ionosphere falls below the frequency of medium-wave 
radio signals allowing them to be transmitted through the ionosphere into outer 
space. The ionosphere appears almost completely transparent to high frequency 
signals such as TV and FM radio signals. Thus, this type of signal is not reflected 





(4.39) 



c 
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by the ionosphere. Consequently, to receive such signals it is necessary to be in 
the line of sight of the relevant transmitter. 



4.6 Faraday rotation 

The electromagnetic force acting on an electron is given by 

f = -e{E + v AB). (4.40) 

If the E and B fields in question are due to an electromagnetic wave propagating 
through a dielectric medium then 

n 

\B\ = -\E\. (4.41) 
c 

It follows that the ratio of the magnetic to the electric forces acting on the elec- 
tron is nv/c. In other words, the magnetic force is completely negligible unless 
the wave amplitude is sufficiently high that the electron moves relativistically in 
response to the wave. This state of affairs is rare, but can occur when intense 
laser beams are made to propagate through plasmas. 

Suppose, however, that the dielectric medium contains an externally gener- 
ated magnetic field B. This can easily be made much stronger than the optical 
magnetic field. In this case, it is possible for a magnetic field to affect the propa- 
gation of low amplitude electromagnetic waves. The electron equation of motion 
(4.11) generalizes to 

ms + fs = -e(E + sAB), (4.42) 

where any damping of the motion has been neglected. Suppose that the direction 
of B is in the positive ^-direction, and that the wave propagates in the same 
direction. With these assumptions the E and s vectors lie in the x-y plane. The 
above equation reduces to 

{ujq^ — uj'^) Sx — iujH Sy = Ex, (4.43a) 

{uq^ - u'^) Sy -\- iujf2 Sx = Ey, (4.43b) 
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provided that all perturbed quantities have an exp(— iwt) time dependence. Here, 

eB 

n = — (4.44) 
m 

is the electron cyclotron frequency. Let 

E± = E^±iEy, (4.45) 

and 

s± = Sx ±isy. (4.46) 

Note that 

= 1{E+ + E_), {4 Ala) 



2 

Ey = 1(E+-E_). (4.47b) 



Equations (4.43) reduce to 

{uq^ - u'^ - u f2) s+ = -—E+, (4.48a) 



m 
e 



(wq - + a; 12) s_ = E_. (4.48b) 



m 



Defining P± = Px±iPy, it follows from Eq. (4.10) that 



(Ne'^/m) E± 



p± = .:;^:r^ - (4-49) 



Finally, from Eq. (4.15), we can write 

e± = n±2 = l + ^, (4.50) 

giving 

2/ ^ -, (iVe^/eom) /. 
n^^cj =1 + — ^ ^ 4.51 
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According to the dispersion relation (4.51), the refractive index of a magne- 
tized dielectric medium can take one of two possible values, which presumably 
correspond to two different types of wave propagating along the ^-axis. The first 
wave has the refractive index n+ and an associated electric field [see Eqs. (4.45)] 

Ex = Eocos[{uj/c){n+z - ct)], (4.52a) 
Ey = Eosm[{uj/c){n+z - ct)]. (4.52b) 

This corresponds to a left-handed circularly polarized wave propagating in the z- 
direction with the phase velocity c/n+. The second wave has the refractive index 
n_ and an associated electric field 

Ex = EQCos[{uj/c){n-z — ct)], (4.53a) 
Ey = -Eosm[{uj/c){n_z - ct)]. (4.53b) 

This corresponds to a right-handed circularly polarized wave propagating in the z- 
direction with the phase velocity c/u-. It is clear from Eq. (4.51) that > ri-. 
Thus, we conclude that in the presence of a ^-directed magnetic field, a ^-directed 
left-handed circularly polarized wave propagates with a phase velocity which is 
slightly less than that of the corresponding right-handed wave. It should be 
remarked that the refractive index is always real (in the absence of damping) , so 
the magnetic field gives rise to no net absorption of electromagnetic radiation. 
This is not surprising since the magnetic field does no work on charged particles, 
and can therefore transfer no energy to the particles from any waves propagating 
through the medium. 

We have seen that right-handed and left-handed circularly polarized waves 
propagate with different phase velocities through a magnetized dielectric medium. 
But, what does this imply for the propagation of a plane polarized wave? Let us 
superimpose the left-handed wave whose electric field is given by Eqs. (4.52) on 
the right-handed wave whose electric field is given by Eqs. (4.53). In the absence 
of a magnetic field n+ = n_ = n, and we obtain 

Ex = 2Eocos[(uj/c){nz - ct)], (4.54a) 
Ey = 0. (4.54b) 
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This, of course, is the field of a plane polarized wave (polarized along the x- 
direction) propagating along the ^-axis with the phase velocity c/n. In the pres- 
ence of a magnetic field we obtain 

Ex = 2£^o cos[(a;/c)(n^ — ct)] cos[(a;/2c)(n+ — n_)^], (4.55a) 

Ey = 2Eqcos[{uj / c){nz — ct)]s\n[{uj/2c){nj^ — n-)z\, (4.55b) 

where 

n=-(n++n_) (4.56) 

is the mean index of refraction. Equations (4.55) can be recognized as the field 
of a plane polarized wave whose angle of polarization with respect to the x-axis, 

X = tan-i(E^/E,), (4.57) 

rotates as the wave propagates along the ^-axis with the phase velocity c/n. In 
fact, the angle of polarization is given by 

X = ^(n+-n_)^, (4.58) 

which clearly increases linearly with the distance traveled by the wave along the 
direction of the magnetic field. This rotation of the plane of polarization of a 
linearly polarized wave propagating through a magnetized dielectric medium is 
known as Faraday rotation (since it was discovered by Michael Faraday in 1845). 

Assuming that the cyclotron frequency Q is relatively small compared to the 
wave frequency cj, and also that uj does not lie close to the resonant frequency 
Wo, it is easily demonstrated that 

(4.59) 

and 

n+-n_~ — ^ 2\2- (^-^0) 

tQmn [uJq — uj'^Y 

It follows that the rate at which the plane of polarization of an electromagnetic 
wave rotates with the distance traveled by the wave is given by 

dx k{u;) NBu 



dl n{uj) 



(4.61) 
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where 5|| is the component of the magnetic field along the direction of propagation 
of the wave, and 



^H= o, ^2....2 .,2>2 - (4-62) 



If the medium in question is a tenuous plasma then n ~ 1 and a;o = 0. Thus, 

dx NB 



dl 2tQm?c o;^ 



(4.63) 



Clearly, the rate at which the plane of polarization rotates is proportional to the 
product of the electron number density and the parallel magnetic field strength. 
Moreover, the plane of rotation rotates faster for low frequency waves than for 
high frequency waves. The total angle by which the plane of polarization is 
twisted after passing through a magnetized plasma is given by 

- o 1 2 / ^(0^11 (0 (4.64) 

provided that N and B\\ vary on length-scales which are large compared to the 
wavelength of the radiation. This formula is regularly employed in radio astron- 
omy to infer the magnetic field-strength in interstellar space. 



4.7 Wave propagation through a magnetized plasma 

For a plasma {ujq = 0) the dispersion relation (4.51) reduces to 

nl(uj) = 1- ^ (4.65) 

The upper sign corresponds to a left-handed circularly polarized wave and the 
lower sign to a right-handed polarized wave. Of course, Eq. (4.65) is only valid 
for wave propagation along the direction of the magnetic field. Wave propagation 
through the Earth's ionosphere is well described by the above dispersion relation. 
There are wide frequency intervals where one of or is positive and the other 
negative. At such frequencies one state of circular polarization cannot propagate 
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through the plasma. Consequently, a wave of that polarization incident on the 
plasma is totally reflected. The other state of polarization is partially transmitted. 

The behaviour of n_^{u) at low frequencies is responsible for a strange phe- 
nomenon known to radio hams as "whistlers." As the frequency tends to zero, 
Eq. (4.65) yields 



2 

n ~ — ^ 



At this sort of frequency is negative, so only right-hand polarized waves can 
propagate. The wave-number of such waves is given by 

Energy transport is governed by the group velocity (see later) 

Thus, low frequency waves transmit energy slower than high frequency waves. 
A lightning strike in one hemisphere of the Earth generates a wide spectrum of 
radiation, some of which propagates along the dipolar field lines of the Earth's 
magnetic field in a manner described approximately by the dispersion relation 
(4.68). The high frequency components of the signal return to the surface of the 
Earth before the low frequency components (since they travel faster along the 
magnetic field) . This gives rise to a radio signal which begins at a high frequency 
and then "whistles" down to lower frequencies. 



4.8 The propagation of electromagnetic radiation through 
a dispersive medium 

Let us now investigate the propagation of electromagnetic radiation through a 
dispersive medium by studying a simple one-dimensional problem. Suppose that 
our dispersive medium extends from x = 0, where it interfaces with a vacuum, 
to a: = oo. Suppose further that a wave is incident normally on the medium. 
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so that the field quantities only depend on x and t. The wave is specified as a 
given function of t at a:; = 0. Since we are not interested in the refiected wave, 
let this function, f{t), say, give the wave amplitude just inside the surface of the 
dispersive medium. Suppose that the wave arrives at this surface at t = 0, and 
that 



How does the wave subsequently develop in the region a; > 0? In order to answer 
this question we must first of all decompose f{t) into harmonic components of 
the form exp(— iwt) {i.e., Fourier harmonics). Unfortunately, if we attempt this 
using only real frequencies, uj, we encounter convergence difficulties, since f{t) 
does not vanish at t = oo. For the moment, we can circumvent these difficulties 
by only considering finite (in time) wave forms. In other words, we now imagine 
that f{t)=0 for t < and t > T. Such a wave form can be thought of as the 
superposition of two infinite (in time) wave forms, the first beginning at t = 
and the second at t = T with the opposite phase, so that the two cancel for all 
time t > T. 

According to standard Fourier transform theory 



Since f{t) is a real function of t which is zero for t < and t > T, we can write 




for t < 0, 
for t > 0. 



(4.69) 



-1 POO POO 

/W = ^ J J J{t')e 



-ia;(t-t')^^/. 



(4.70) 




(4.71) 




(4.72) 



Equation (4.69) yields 



■j^ POO pi / 27rt' \ 
f{t) = - duo I sin I ) cos[u{t - t') ] dt' , 



Jo Jo 




(4.73) 
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or 

_ J_ f°° cos[27rt7r + ujjt - t') ] _ cos[27rt7T - uj{t - t') ] 
~ 2^ Jo ^\ uj - 27r/T uj + 27r/T 

(4.74) 

Let us assume, for the sake of simplicity, that 

T = Nt, (4.75) 

where is a positive integer. This ensures that f{t) is continuous at t = T. 
Equation (4.74) reduces to 

2 /'°° duj 

^^^^ ^tJo ^2 _ (2^/^)2 - T)] - cosu;0 . (4.76) 

This expression can be written 




2 (cos[a;(t-T)] -coswt), (4.77) 



or 

It is not entirely obvious that Eq. (4.78) is equivalent to Eq. (4.77). However, we 
can easily prove that this is the case by taking Eq. (4.78) and using the standard 
definition of a real part (i.e., half the sum of the quantity in question and its 
complex conjugate) to give 

- i /1 -'i 

Replacing the dummy integration variable u by —u in the second integral and 
then making use of symmetry, it is easily seen that the above expression reduces 
to Eq. (4.77). 
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CO = 27C/T 



Figure 6: Sketch of the integration contours used to evaluate Eqs. (4.78) and 
(4.81) 

Equation (4.77) can be written 

m = 1 11 ^ sin[.(t - T/2)] (4.80) 

Note that the integrand is finite at a; = 27r/r, since at this point the vanishing of 
the denominator is compensated for by the simultaneous vanishing of the numer- 
ator. It follows that the integrand in Eq. (4.78) is also not infinite at a; = 27r/r, 
as long as we do not separate the two exponentials. Thus, we can replace the 
integration along the real axis through this point by a small semi-circle in the 
upper half of the complex plane. Once this has been done, we can deform the 
path still further and can integrate the two exponentials in Eq. (4.78) separately: 

fit) = ^Ref e-^-^ t / - T^Re / e-^-('-^^—^ (4.81) 
^ ^ 27r Jc uj- 27r/r 27t Jc uj - 27r/r ^ ^ 

The contour C is sketched in Fig. 6. Note that it runs from +oo to — oo, which 
accounts for the change of sign between Eqs. (4.78) and (4.81). 

We have already noted that a finite wave form which is zero for t < and 
t > T can be through of as the superposition of two out of phase infinite wave 
forms, one starting at t = and the other at t = T. It is plausible, therefore, 
that the first term in the above expression corresponds to the infinite wave form 
starting at t = 0, and the second to the infinite wave form starting at t = T. If 
this is the case then the signal (4.69), which starts at t = and ends at t = oo. 



147 



can be written in the form 



(4.82) 



a; — 27r/r 



Let us test this proposition. In order to do this we must replace the original path 
of integration C by two equivalent paths. 

First, consider t < 0. In this case, —\ujt has a negative real part in the upper 
half plane which increases indefinitely with increasing distance from the axis. 
Thus, we can replace the original path of integration by the path A (see Fig. 7). 
The integral clearly vanishes along this path if we let A approach infinity in the 
upper half plane. Consequently, 



Next, consider t > 0. Now, —loot has a negative real part in the lower half 
plane, so that the exponential vanishes at infinity in this half plane. If we attempt 
to deform C to infinity in the lower half plane, the path of integration "catches" on 
the singularity of the integrand at a; = 27r/r (see Fig. 7). The path of integration 
B therefore consists of three parts: the part at infinity, Bi, where the integral 
vanishes due to the exponential factor e~"^*; B2, the two parts leading to infinity 
which cancel each other and thus contribute nothing to the integral; the path ^3 
around the singularity. This latter contribution can easily be evaluated using the 
Cauchy residue theorem: 



Thus, it is proven that the expression (4.82) actually describes a wave form be- 
ginning at t = whose subsequent motion is specified by Eq. (4.69). 

Equation (4.82) can immediately be generalized to give the wave motion in 
the region a; > 0: 



fit) = 



(4.83) 



for t < 0. 




(4.84) 




(4.85) 
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A 




Figure 7: Sketch of the integration contours used to evaluate Eq. (4.82) 

This follows from standard wave theory, because we know that an unterminated 
wave motion at a; = of the form e~"^* takes the form ^'^i^^-^^) after moving a 
distance x in the dispersive medium, provided that k and uj are related by the 
appropriate dispersion relation. For a medium consisting of a single resonant 
species this dispersion relation is written (see Eq. (4.17) ) 




uS^ — igujujo 



(4.86) 



4.9 Propagation of the wave front in a dispersive medium 

It is helpful to define 

s = t--. (4.87) 

c 

Let us consider the two cases s < and s > separately. 

Suppose that s < 0. In this case we distort the path C, used to evaluate the 
integral (4.85), into the path A shown in Fig. 8. This is only a sensible thing to 
do if the real part ofi{kx — cot) is negative at infinity in the upper half plane. It 
is clear from the dispersion relation (4.86) that k = uj/c in the limit — > oo. 
Thus, 

i {kx — ujt) = —'\uj{t — x/c) = —iujs. (4.88) 
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It follows that i [kx—ujt) possesses a large negative real part along path A provided 
that s < 0. Thus, Eq. (4.85) yields 



f(x,t) = {) (4.89) 

for s < 0. In other words, it is impossible for the wave front to propagate through 
the dispersive medium with a velocity greater than the velocity of light in a vacuum. 

Suppose that s > 0. In this case we distort the path C into the lower half 
plane, since i {kx—cot) = — i cos has a negative real part at infinity in this region. In 
doing this, the path becomes stuck not only at the singularity of the denominator 
when UJ = 2'k/t, but also at the branch points of the expression for k. After a 
little algebra, the dispersion relation (4.86) yields 

(4.90) 

C y UJQ-\. — UJ \j UJQ— — UJ 

where 

a;o± = -ip± ^JuQ^ -p^, (4.91) 

and 




= -i p ± ^uj^^ + a;/ - p2 . (4.92) 

Here, 

^Ne'^/tom (4.93) 



UJp 



is the plasma frequency, and 



p=3^^^, (4.94) 



parameterizes the damping. In order to prevent multiple roots of Eq. (4.90) it is 
necessary to place branch cuts between a;o+ and and also between ujq_ and 
uj\- (see Fig. 8). 

The path of integration B is conveniently split into the parts Bi through 
B^. The contribution from Bi is negligible since the exponential in Eq. (4.85) is 
vanishingly small on this part of the integration path. Likewise, the contribution 
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(0 = 271/X. 







branch cu\ 





Figure 8: Sketch of the integration contours used to evaluate Eq. (4.85) 



from B2 is zero since its two sections always cancel. The contribution from ^3 
follows from the residue theorem: 



^3 = J-Re(27rie^['=-^-2"*/^]). 



(4.95) 



Here, kr denotes the value of k obtained from the dispersion relation (4.86) in 
the limit uj — > 27: /t. Thus, 



53 = e-^^^C^-)^ sin (^t: ^ - Re(kr) 



(4.96) 



In general, the contributions from B4 and cannot be simplified further. For 
the moment we denote them as 



Ba = — Re i eHkx-u;t), 



and 



27r T^^ UJ — 2'k/t^ 



27r Jb uj - 27r/r 



(4.97) 



(4.98) 
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where the paths of integration circle the appropriate branch cuts. In all, we have 
f{x, t) = e-i^^C^-) ^ sin ^27r ^ - Re(fc^) + B4 + B5 (4.99) 

for s > 0. 

Let us now look at the special case s = 0. For this value of s we can change the 
original path of integration to one at infinity in either the upper or the lower half 
plane, since the integrand vanishes in each case, through no longer exponentially, 
but rather as l/w^. We can see this from Eq. (4.82), which can be written in the 
form 

fit) = ^(f e-'"'^i- + / e+'-'^V) • (4-100) 
Substitution of uj for —uj in the second integral yields 

Now, applying dispersion theory, we get from the above equation, just as we got 
Eq. (4.85) from Eq. (4.82), 

/(.,*) = i/e'(—)^^,-|^. (4.102) 

Clearly, the integrand vanishes as e~^'^^/a;^ as uj becomes very large. Thus, it 
vanishes as l/o;^ for s = 0. Since we can calculate f{x,t) by using either path A 
or path B, we can see that 



f{x,t) = e-^'"^'^-^)^ sin^27r ^ - Re(A;^) x j + ^4 + ^5 = 



(4.103) 



for s = 0. Thus, there is continuity in the transition from the region s < to the 
region s > 0. 

We are now in a position to make some meaningful statements about the 
behaviour of the signal at depth x inside the dispersive medium. Prior to the 
time t = x/c there is no motion. Even if the phase velocity is superluminal, no 
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electromagnetic signal can arrive earlier than one propagating with the velocity of 
light in vacuum c. The wave motion for t > x/c is conveniently divided into two 
parts: free oscillations and forced oscillations. The former are given by 54 + ^5, 
and the latter by 



X 

t 



(4.104) 



where 



is termed the phase velocity. The forced oscillations have the same sine wave 
characteristics and oscillation frequency as the incident wave. However, the wave 
amplitude is diminished by the damping coefficient, although, as we have seen, 
this is generally a negligible effect unless the frequency of the incident wave closely 
matches one of the resonant frequencies of the dispersive medium. The phase ve- 
locity Vp determines the velocity with which a point of constant phase {e.g., a peak 
or trough) of the forced oscillation signal propagates into the medium. However, 
the phase velocity has no effect on the velocity with which the forced oscillation 
wave front propagates into the medium. This latter velocity is equivalent to the 
velocity of light in vacuum c . The phase velocity Vp can be either greater or 
less than c, in which case peaks and troughs either catch up with or fall further 
behind the wave front. Of course, peaks can never overtake the wave front. 

It is clear from Eqs. (4.91), (4.92), (4.97), and (4.98) that the free oscilla- 
tions oscillate with real frequencies which are somewhere between the resonant 
frequency luq and the plasma frequency ujp. Furthermore, the free oscillations are 
damped in time like ex.p{—pt). The free oscillations, like the forced oscillations, 
begin at time t = x/c. At t = x/c the free and forced oscillations just cancel (see 
Eq. (4.103) ). As t increases both the free and forced oscillations set in, but the 
former rapidly damp away, leaving only the forced oscillations. Thus, the free 
oscillations can be regarded as some sort of transient response of the medium to 
the incident wave, whereas the forced oscillations determine the time asymptotic 
response. The real frequency of the forced oscillations is that imposed exter- 
nally by the incident wave, whereas the real frequency of the free oscillations is 
determined by the nature of the dispersive medium, quite independently of the 
frequency of the incident wave. 
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One slightly surprising result of the above analysis is the prediction that the 
wave front of the signal propagates into the dispersive medium with the veloc- 
ity of light in vacuum, irrespective of the dispersive properties of the medium. 
Actually, this is a fairly obvious result. As is well described by Feynman in his 
famous Lectures on Physics, when an electromagnetic wave propagates through a 
dispersive medium, the electrons and ions which make up that medium oscillate 
in sympathy with the incident wave and in doing so emit radiation. Both the 
radiation from the electrons and ions and the incident radiation travel at the 
velocity c. However, when these two radiation signals are superposed the net 
effect is as if the incident signal propagates through the dispersive medium with 
a phase velocity which is different from c. Consider the wave front of the incident 
signal, which clearly propagates into the medium with the velocity c. Prior to the 
arrival of this wave front the electrons and ions are at rest, since no information 
regarding the arrival of the incident wave at the surface of medium can propagate 
faster than c. After the arrival of the wave front the electrons and ions are set 
into motion and emit radiation which can affect the apparent phase velocity of 
radiation which arrives somewhat later. But this radiation certainly cannot affect 
the propagation velocity of the wave front itself, which has already passed by the 
time the electrons and ions are set into motion (because of the finite inertia of 
the electrons and ions). 

4.10 The Sommerfeld precursor 

Let us consider the situation immediately after the arrival of the signal; i.e., when 
s is small and positive. Let us start from Eq. (4.102), which can be written in 
the form 

= i/^e'«-Wc,.-..)__^. (4.106) 

We can deform the original path of integration C into a large semi-circle of radius 
R in the upper half-plane, plus the segments of the real axis, as shown in Fig. 9. 
Because of the denominator — (27r/T)^, the integrand tends to zero as 1/a;^ 
on the real axis. We may add the path in the lower half plane which is shown 
as a dotted line in the figure, for if the radius of the semi-circular portion of this 
lower path is increased to infinity, the integrand vanishes exponentially because 
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Figure 9: Sketch of the integration contour used to evaluate Eq. (4.107) 

s > 0. Therefore, we may replace our original path of integration by the entire 
circle S. Thus, 

= 1 f/'-"'-'''- \._%,^y (4.107) 
in the limit that the radius of the circle R tends to infinity. 
The dispersion relation (4.86) yields 



p 



up' I 2 c a; 



(4.108) 



in the limit — > oo. Using the abbreviation 

2 



P 



2c 



(4.109) 



and henceforth neglecting 27: /r with respect to w, we obtain from Eq. (4.107) 



f{x,t) = fi(^,t) ~ i / exp 

T Js 



duj 



(4.110) 
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This expression can also be written 



fi(C,t) = i / exp 
T Js 



— 1 



^. (4.111) 



UJ 



Let 



1 u 



It follows that 



giving 



— = 1 du, 

UJ 



duj Is 



-'""du. 



(4.112) 



(4.113) 



(4.114) 



Substituting the angular variable u for uj as the integration variable in Eq. (4.111) 
yields 



1 s 



h{i,t) = - J- j exp(-2iV^s cosw)e '"du. 



(4.115) 



Here, we have taken \/^/s as the radius of the circular integration path in the 
w-plane. This is indeed a large radius, since s <C 1. From symmetry, Eq. (4.115) 
simplifies to 

f\{itt) = — \ -z / exp(— 2iy^s cosu) cosudu. (4.116) 
V wo 



The following mathematical identity is very well-known^ ^ 



II f-1-K 



A z cos Q 



27r 



cos(n^) dQ^ 



(4.117) 



where Jn(-z) is Bessel function of order n. It follows from Eq. (4.115) that 



/iK,t) = ^y|ji(2^). 



(4.118) 
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Figure 10: The Bessel function Ji{z) 



Here, we have made use of the fact that Ji{—z) = —Ji{z). 

The properties of Bessel functions are well-known and are listed in many 
standard references on mathematical functions (see, for instance, Abramowitz 
and Stegun). In the small argument limit ^ <C 1 we find that 



The behaviour of Ji{z) is further illustrated in Fig. 10. 

We are now in a position to make some quantitative statements regarding 
the signal which first arrives at depth x in the dispersive medium. This sig- 
nal propagates at the velocity of light in vacuum and is called the Sommerfeld 

^^M. Abramowitz, and LA. Stegun, Handbook of mathematical functions, (Dover, New York, 
1965), Eq. 9.1.21. 



Mz) = '- + Oiz% 



(4.119) 



On the other hand, in the large argument limit ^ 1 we obtain 




(4.120) 
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Figure 11: The Sommerfeld precursor 

precursor. The first important point to note is that the amphtude of the Som- 
merfeld precursor is very small compared to that of the incident wave (whose 
amplitude is normalized to unity). We can easily see this because in deriving 
Eq. (4.118) we assumed that |a;| = y^^/s 27r/r on the circular integration 
path S. Since the magnitude of Ji is always less than, or of order, unity, it is 
clear that <C 1. This is a comforting result, since in a naive treatment of 
wave propagation through a dielectric medium the wave front propagates at the 
group velocity Vg (which is usually less than c) and, therefore, no signal should 
reach depth x in the medium before time x/vg. We are finding that there is, in 
fact, a precursor which arrives at t = x/c, but that this signal is fairly small. 
Note from Eq. (4.109) that ^ is proportional to x. Clearly, the amplitude of the 
Sommerfeld precursor decreases like one over the distance traveled by the wave 
front through the dispersive medium (since Ji attains its maximum value when 
s ~ 1/0- Thus, the Sommerfeld precursor is likely to become undetectable after 
the wave has traveled a long distance through the medium. 




158 



Equation (4.118) can be written 



fi{U) = ^9{s/so), (4.121) 

where sq = 1/4^, and 

g{z) = V^Ji{V^). (4.122) 

The normahzed Sommerfeld precursor g[z) is shown in Fig. 11. It can be seen that 
both the ampHtude and the oscillation period of the precursor gradually increase. 
The roots of Ji{z) [i.e., the solutions of Ji{z) = 0] are spaced at distances of 
approximately tt apart. Thus, the time interval for the mth half period of the 
precursor is approximately given by 

2 

777 TT 

Atm ~ (4.123) 
Note that the initial period of oscillation, 

Ato ~ ^, (4.124) 

is extremely small compared to the incident period r. Moreover, the initial period 
of oscillation is completely independent of the frequency of the incident wave. In 
fact, Ato depends only on the depth x and on the dispersive power of the medium. 
The period decreases with increasing distance x traveled by the wave front though 
the medium. So, when visible radiation is incident on some dispersive medium it 
is quite possible for the first signal detected well inside the medium to lie in the 
X-ray region of the electromagnetic spectrum. 



4.11 The method of stationary phase 

Equation (4.102) can be written in the form 

f{x,t)= f e''^^^^F{uj)duj (4.125) 
Jc 



159 



where 

^H = -^— rW^' (4-126) 

and 

(t){uj) = k{uj)x -ujt. (4.127) 

It is clear that F[uj) is a relatively slowly varying function of uo (except in the 
immediate vicinity of the singular points uj = ±27r/r), whereas the phase (j){uo) is 
generally large and rapidly varying. The rapid oscillations of exp(i0) over most 
of the range of integration means that the integrand averages to almost zero. 
Exceptions to this cancellation rule occur only when (j){uS) is stationary] i.e., when 
4>{uS) has an extremum. The integral can therefore be estimated by finding places 
where 4>{uS) has a vanishing derivative, evaluating (approximately) the integral in 
the neighbourhood of each of these points, and summing the contributions. This 
procedure is called the method of stationary phase. 

Suppose that 4'{oo) has a vanishing first derivative at a; = a;^. In the neigh- 
bourhood of this point, (j){uS) can be expanded as a Taylor series, 

0(a;) = 0, + i0';(u; - ujsf + • • • . (4.128) 

Here, the subscript s is used to indicate or its second derivative evaluated at 
UJ = uOg- Since F{uj) is slowly varying, the contribution to the integral from this 
stationary phase point is approximately 



-oo 



f,c^F{us)e"^^ / Qm)<t>s{^-^s) dcj. (4.129) 

J oo 

It is tacitly assumed that the stationary point lies on the real axis in a;-space, so 
that locally the integral along the contour C is an integral along the real axis in 
the direction of decreasing cu. The above expression can be written in the form 

I A POO 

fs ~ -F(a;,) e'^^M ^ [cos(7rtV2) + i sin(7rtV2)] dt, (4.130) 



where 



lt' = l,l>':{u>-u;,)'. (4.131) 
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The integrals in the above expression are, Fresnel integrals^^ and can be shown 
to take the values 

/>oo />00 -1 

/ cos(7rtV2) dt= sin(7rtV2) dt= -. (4.132) 
Jo Jo 2 

It follows that 



Ac.-J^FK)e'*-. (4.133) 

It is easily seen that the arc length (in cj-space) of the integration contour which 
makes a significant contribution to fg is of order Acu/cUs ~ 1/ ^Jk{ujs) x. Thus, 
the arc length is relatively short provided that the wavelength of the signal is 
much less than the distance propagated through the dispersive medium. If there 
is more than one point of stationary phase in the range of integration then the 
integral is approximated as a sum of terms like the above. 

Integrals of the form (4.125) can be calculated exactly using the method of 
steepest decent}^ The stationary phase approximation (4.133) agrees with the 
leading term of the method of steepest decent (which is far more difficult to 
implement than the method of stationary phase) provided that 4>{uS) is real [i.e., 
provided that the stationary point lies on the real axis). If is complex, however, 
the stationary phase method can yield erroneous results. This suggests that the 
stationary phase method is likely to break down when the extremum point uj = LOg 
approaches any poles or branch cuts in the w-plane (see Fig. 8). 



4.12 The group velocity 

The point of stationary phase, defined by d4>/duj = 0, satisfies the condition 

c ct , . 

- = -, 4.134 

v„ X 



^^M. Abramowitz, and LA. Stegun, Handbook of mathematical functions, (Dover, New York, 
1965), Sec. 7.3. 

^^Leon Brillouin, Wave propagation and group velocity, (Academic press. New York, 1960). 
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where 



V, 



duj 
dk 



(4.135) 



is conventionally termed the group velocity. Thus, the signal seen at position x 
and time t is dominated by the frequency range whose group velocity Vg is equal 
to x/t. In this respect, the signal incident at the surface of the medium [x = 0) 
at time t = can be said to propagate through the medium at the group velocity 

Vg{uj). 



The simple one-resonance dispersion relation (4.86) yields 



— ~ n{Lxj) 



V, 



1 + 



+ 



UJ 



UJr 



UJ' 



UJ^ — UJr 



UJ. 



P J 



(4.136) 



in the limit g 0, where 



n{uj] = — 

UJ 



P 



UJ' 



UJq — Uj'^ 



(4.137) 



The variation of c/vg and the refractive index n with frequency is sketched in 
Fig. 12. With ^ = the group velocity is less than c for all a;, except for 

a;o < ^ < ^1 = uJq' + uj^ , where it is purely imaginary. Note that the refractive 

index is also complex in this frequency range. The phase velocity Vp = c/n is 
subluminal for uj <ujq, imaginary for ujq < uj < uj\, and superluminal for uj > uj\. 

The frequency range which contributes to the amplitude at time t is deter- 
mined graphically by finding the intersection of a horizontal line with ordinate 
ct/x with the solid curve in Fig. 12. There is no crossing of the two curves for 
t < to = x/c, thus no signal arrives before this time. For times immediately 
following to the point of stationary phase is seen to be at a; — > oo. In this large 
UJ limit the point of stationary phase is given by 



UJg ~ UJr 



to 



(4.138) 



Note that uj = 
that 



'pit-toY 

-ujg is also a point of stationary phase. It is easily demonstrated 



0, ~ -2^^(t-to): 



(4.139) 
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Figure 12: The typical variation of the functions c/vg{uS) and n{uS). Here, uj\ = 
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and 



with 



F{ujs) - 



(4.140) 



(4.141) 



Here, ^ is given by Eq. (4.109). The stationary phase approximation (4.133) gives 



+ C.C., 



(4.142) 



where c.c. denotes the complex conjugate of the preceding term (this contribution 
comes from the second point of stationary phase located at a; = —ujs ) . The above 
expression reduces to 



f^-— ^57^^°^ 



^V^it - to) - 37r/4 



(4.143) 



It is easily demonstrated that the above formula is the same as the expression 
(4.118) for the Sommerfeld precursor in the large argument limit t — Iq ^ 1/^. 
Thus, the method of stationary phase yields an expression for the Sommerfeld 
precursor which is accurate at all times except those immediately following the 
first arrival of the signal. 



4.13 The Brillouin precursor 

As time progresses the horizontal line ct/x in Fig. 12 gradually rises and the 
point of stationary phase moves to ever lower frequencies. In general, however, 
the amplitude remains relatively small. Only when the elapsed time reaches 

h = > to 4.144 

c 

is there a qualitative change. This time marks the arrival of a second precursor 
known as the Brillouin precursor. The reason for the qualitative change is ev- 
ident from Fig. 12. At t = ti the lower region of the c/vg curve is intersected 
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for the first time, and uj = becomes a point of stationary phase. It is clear 
that the oscillation frequency of the Brillouin precursor is far less than that of 
the Sommerfeld precursor. Moreover, it is easily demonstrated that the second 
derivative of k{uj) vanishes at a; = 0. This means that 0" = 0. The stationary 
phase result (4.133) gives an infinite answer in such circumstances. Of course, the 
amplitude of the Brillouin precursor is not infinite, but it is significantly larger 
than that of the Sommerfeld precursor. 

In order to generalize the result (4.133) to deal with a stationary phase point 
at cj = it is necessary to expand 0(a;) about this point, keeping terms up to cu^. 
Thus, 

^ (4.145) 



0(a;)~a;(^l-^) + ^fc^"a;^ 



where 



k'" ^ 



3uj. 



p 



for the simple dispersion relation (4.86). The amplitude (4.125) is therefore given 
approximately by 

This expression reduces to 



-oo 



(4.147) 



oo 










COS 













dv, 



(4.148) 



where 



V 



and 



z = 




3vW 



(4.149) 



(4.150) 



The positive (negative) sign in the integrand is taken for t <t\ {t > ti). 
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Figure 13: A sketch of the behaviour of the Brillouin precursor as a function of 
time 



The integral in Eq. (4.150) is known as an Airy integral. It can be expressed 
in terms of Bessel functions of order 1/3, as follows: 

,3 N 



-oo 



COS 



6 

{ — + v 



dv = 



and 



poo 
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TT 



d'^= ^ [Jl/3{z) + J-l/3{z)] . 



(4.151) 



(4.152) 



From the well-known properties of Bessel functions the precursor can be seen 
to have a growing exponential character for times earlier than t = ti, and an 
oscillating character for t > ti. The amplitude in the neighbourhood of t = ti is 
plotted in Fig. 13. 

The initial oscillation period of the Brillouin precursor is crudely estimated 
(from ^ ~ 1) as 

Ato ~ {xk'^'y/^. (4.153) 
The amplitude of the Brillouin precursor is approximately 

T 



I/I 



(4.154) 
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Let us adopt the ordering 



l/r ~ Wo ~ <C (4.155) 

which corresponds to most physical situations involving the propagation of elec- 
tromagnetic radiation through dielectric media. It follows from the above results, 
plus the results of Section 4.10, that 

(Atoa;p)briiiouin ~ — > 1, (4.156) 



and 



Furthermore, 



{Ato UJp) sommevfeld ^ \ ^ (4.157) 



/ \ 1/3 

l/lbrillouin ~ ( ^ < 1, (4.158) 



and 



|/|sommerfeld ~ ("^^ ^ |/|brillouin- (4.159) 

It is clear that the Sommerfeld precursor is a low amplitude, high frequency signal, 
whereas the Brillouin precursor is a higher amplitude, low frequency signal. Note 
that the amplitude of the Brillouin precursor, whilst it is significantly higher than 
that of the Sommerfeld precursor, is still much less than that of the incident wave. 



4.14 Signal arrival 

Let us try to establish at what time t2 a signal first arrives at position x inside the 
dielectric medium whose amplitude is comparable with that of the wave incident 
at time t = on the surface of the medium (x = 0). Let us term this event 
the "arrival" of the signal. It is plausible from the discussion in Section 4.11 
regarding the stationary phase approximation that signal arrival corresponds to 
the situation where the point of stationary phase in cj-space corresponds to a 
pole of the function F{uj). In other words, when cus approaches the frequency 
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t -> 



precursor 

Figure 14: A sketch of the signal ampHtude as a function of time as seen inside 
some dielectric medium subject to an incident wave which starts at some specific 
time 



27r/T of the incident signal. It is certainly the case that the stationary phase 
approximation yields a particularly large amplitude signal when LOg — > 27: Jr. 
Unfortunately, as has already been discussed, the method of stationary phase 
becomes inaccurate under these circumstances. However, calculations involving 
the more robust method of steepest decent^^ confirm that in most cases the signal 
amplitude first becomes significant when LOg = 27: Jr. Thus, the signal arrival time 
is 

t2 = , y (4.160) 

where i;^(27r/r) is the group velocity calculated using the frequency of the incident 
signal. It is clear from Fig. 12 that 

to<ti< t2. (4.161) 

Thus, the main signal arrives later than the Sommerfeld and Brillouin precursors. 

The final picture which emerges from our investigations is summarized in 
Fig. 14. The main signal arrives at the group velocity corresponding to the fre- 
quency of the incident wave. However, it is possible to detect the arrival of the sig- 
nal before this, given sufficiently accurate detection equipment. In fact, the first 
information regarding the arrival of the incident wave at the vacuum/dielectric 
interface propagates at the velocity of light in a vacuum. 



14 



Leon Brillouin, Wave propagation and group velocity, (Academic press, New York, 1960). 
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4.15 The propagation of radio waves through the iono- 
sphere 



We have studied the transient behaviour of an electromagnetic wave incident 
on a spatially uniform dielectric medium in great detail. Let us now consider 
a quite different, but equally important, problem. What is the time asymptotic 
steady-state behaviour of an electromagnetic wave propagating though a spatially 
non-uniform dielectric medium? 

As a specific example, let us consider the propagation of radio waves through 
the Earth's ionosphere. The refractive index of the ionosphere can be written [see 
Eq. (4.27)] 

2 -, 



= 1 - , ^ . , , (4.162) 

where is a real positive constant which parameterizes the damping of electron 
motion (in fact, v is the collision frequency of free electrons with other particles 
in the ionosphere), and 

is the plasma frequency. In the above formula, N is the density of free electrons in 
the ionosphere and m is the electron mass. We shall assume that the ionosphere 
is horizontally stratified, so that N = N{z), where the coordinate z measures 
height above the Earth's surface {n.b., the curvature of the Earth is neglected in 
the following analysis). The ionosphere actually consists of two main layers; the E- 
layer, and the F-layer. We shall concentrate on the lower E-layer, which lies about 
100 km above the surface of the Earth, and is about 50 km thick. The typical 
day-time number density of free electrons in the E-layer is 3 x 10^^ m~^. 

At night-time, the density of free electrons falls to about half this number. The 
typical day-time plasma frequency of the E-layer is, therefore, about 5 MHz. The 
typical collision frequency of free electrons in the E-layer is about 0.05 MHz. 
According to simplistic theory, any radio wave whose frequency lies below the 
day-time plasma frequency, 5 MHz, [i.e., any wave whose wavelength exceeds 
about 60 m) is reflected by the ionosphere during the day. Let us investigate 
in more detail exactly how this process takes place. Note, incidentally, that for 
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mega-Hertz frequency radio waves <C so it follows from Eq. (4.162) that 
is predominately real {i.e., under most circumstances, the electron collisions can 
be neglected). 

The problem of radio wave propagation through the ionosphere was of great 
practical importance during the first half of the 20th Century, since at that 
time long- wave radio waves were the principle means of military communication. 
Nowadays, the military have far more reliable ways of communicating. Neverthe- 
less, this subject area is still worth studying because the principle tool used to 
deal with the problem of wave propagation through a non-uniform medium, the 
so-called W.K.B. approximation, is of great theoretical importance. In particular, 
the W.K.B. approximation is very widely used in quantum mechanics (in fact, 
there is a great similarity between the problem of wave propagation through a 
non-uniform medium and the problem of solving Schrodinger's equation in the 
presence of a non-uniform potential). 

Maxwell's equations for a wave propagating through a non-uniform, unmag- 
netized, dielectric medium are: 

V-^; = 0, (4.164a) 
V cB = 0, (4.164b) 
VA^; = ikcB, (4.164c) 
VAcB = -ikn'^E, (4.164d) 

where n is the non-uniform refractive index of the medium. It is assumed that all 
field quantities vary in time like e~"^*, where uj = kc. Note that, in the following, 
k is the wavenumber in free space, rather than the wavenumber in the dielectric 
medium. 



4.16 The W.K.B. approximation 

Consider a radio wave which is vertically incident, from below, on the horizontally 
stratified ionosphere. Since the wave normal is initially aligned along the ^-axis, 
and since n = n(^), we expect all field components to be functions of z only, so 
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that 

In this situation, Eqs. (4.164) reduce to = cB^ = 0, with 



dz 
dcBx 
dz 

and 

dE, 



ikcB^, (4.166a) 
-ikn'^Ey, (4.166b) 



dz 
dcBy 
dz 



ikcBy, (4.167a) 
-ikn^E^. (4.167b) 



Note that Eqs. (4.166) and (4.167) are isomorphic and completely independent 
of one another. It follows that, without loss of generality, we can assume that the 
wave is linearly polarized with its electric vector parallel to the y-axis. This means 
that we are only going to consider the solution of Eqs. (4.166). The solution of 
Eqs. (4.167) is of exactly the same form, except that it describes a linear polarized 
wave with its electric vector parallel to the 

Equations (4.166) can be combined to give 

(PE 



^k^T?Ey = 0. (4.168) 



dz'' 



Since Ey is a function of z only, we now use the total derivative sign d/dz instead 
of the partial derivative sign d/dz. The solution of the above equation for the 
case of a uniform medium, where n is constant, is straightforward: 

Ey = Ae''^^'\ (4.169) 

where ^4 is a constant, and 

(t) = ±knz. (4.170) 
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Note that the e~"^* time dependence of all wave quantities is taken as read 
during this investigation. The solution (4.169) represents a wave of constant 
amplitude A and phase (j){z). According to Eq. (4.170), there are, in fact, two 
independent waves which can propagate through the medium in question. The 
upper sign corresponds to a wave which propagates vertically upwards, and the 
lower sign corresponds to a wave which propagates vertically downwards. Both 
waves propagate with the constant phase velocity c/n. 

In general, if n = n{z) the solution of Eq. (4.168) does not remotely resemble 
the wave-like solution (4.169). However, in the limit in which n{z) is a "slowly 
varying" function of z (exactly how slowly varying is something which we shall 
establish later), we expect to recover wave-like solutions. Let us suppose that 
n{z) is indeed a "slowly varying" function, and let us try substituting the wave 
solution (4.169) into Eq. (4.168). We obtain 

This is a non-linear differential equation which, in general, is very difficult to 
solve. However, we note that if n is a constant then (P(f)/dz'^ = 0. It is, therefore, 
reasonable to suppose that if n{z) is a "slowly varying" function then the last 
term on the right-hand side of the above equation can be regarded as being small. 
Thus, to a first approximation Eq. (4.171) yields 

^~±A;n, (4.172) 
dz 

and 

It is clear from a comparison of Eqs. (4.171) and (4.173) that n{z) can be regarded 
as a "slowly varying" function of z as long as its variation length-scale is far longer 
than the wavelength of the wave. In other words, provided that (dn/dz) / {k n^) <C 
1. 

The second approximation to the solution is obtained by substituting Eq. (4.173) 
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into the right-hand side of Eq. (4.171): 



This gives 



dz \ dz J 



d(f) , , / , i dnV^^ , , dn 
— ^±kn[l±—^—\ c^±kn+— — , 
dz \ kn-^ dz J In dz 



1/2 



(4.174) 



(4.175) 



where use has been made of the binomial expansion. The above expression can 
be integrated to give 



~ ibfc j ndz + \\og{p}^^). 
Substitution of Eq. (4.176) into Eq. (4.169) yields the final result 

^,.^n-V^exp(±i./n..). 

It follows from Eq. (4.166a) that 



(4.176) 



(4.177) 



cBx ~ ^An^^^ exp i±ik ndz 



i A dn 



2fcn3/2 dz 



— exp 



^±i/c j ndz^ . 



(4.178) 



Note that the second term is small compared to the first, and can usually be 
neglected. 

Let us test to what extent the expression (4.177) is a good solution of Eq. (4.168) 
by substituting this expression into the left-hand side of the equation. The result 
is 

This must be small compared with either term on the left-hand side of Eq. (4.168). 
Hence, the condition for Eq. (4.177) to be a good solution of Eq. (4.168) becomes 



1 



3 f 1 dn\ 



1 d'^n 



4 dz J 2rfi dz^ 



< 1. 



(4.180) 
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The solutions 



. ^.-V^exp(±i./^..), (4.181a) 

cB. . ^^n^/^exp(±i./n..), (4.181b) 

to the non-uniform wave equations (4.166) are most commonly called the W.K.B. 
solutions, in honor of G. Wentzel, H.A. Kramers, and L. Brillouin, who are cred- 
ited with independently discovering these solutions (in a quantum mechanical 
context) in 1926. Actually, H. Jeffries wrote a paper on these solutions (in a 
wave propagation context) in 1923. Hence, some people call these the W.K.B. J. 
solutions (or even the J. W.K.B. solutions). In fact, these solutions were first dis- 
cussed by Liouville and Green in 1837, and again by Rayleigh in 1912. We shall 
refer to Eqs. (4.181) as the W.K.B. solutions, since this is what they are most 
commonly called. However, it should be understand that, in doing so, we are 
not making any statement as to the credit due to various scientists in discovering 
these solutions. After all, this is not a history of science course! 

Recall, that when a propagating wave is normally incident on an interface, 
where the refractive index suddenly changes (for instance, when a light wave 
propagating in the air is normally incident on a glass slab), there is generally 
significant reflection of the wave. However, according to the W.K.B. solutions 
(4.181), when a propagating wave is normally incident on a medium in which the 
refractive index changes slowly along the direction of propagation of the wave, 
then the wave is not reflected at all. This is true even if the refractive index varies 
very substantially along the path of propagation of the wave, as long as it varies 
slowly. The W.K.B. solutions imply that as the wave propagates through the 
medium its wavelength gradually changes. In fact, the wavelength at position 
z is approximately X{z) = 27T/kn{z). Equations (4.181) also imply that the 
amplitude of the wave gradually changes as it propagates. In fact, the amplitude 
of the electric field component is inversely proportional to n^/^, whereas the 
amplitude of the magnetic field component is directly proportional to n^/^. Note, 
however, that the energy flux in the ^-direction, given by the the Poynting vector 
— {EyB^* + Ey* B^)/ {4:/j,q), remains constant (assuming that n is predominately 
real). 
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Of course, the W.K.B. solutions (4.181) are only approximations. In reality, a 
wave propagating into a medium in which the refractive index is a slowly varying 
function of position is subject to a small amount of reflection. However, it is easily 
demonstrated that the ratio of the reflected amplitude to the incident amplitude 
is of order [dn / dz) / {k n'^) . Thus, as long as the refractive index varies on a much 
longer length-scale than the wavelength of the radiation, the reflected wave is 
negligibly small. This conclusion remains valid as long as the inequality (4.180) 
is satisfied. There are two main reasons why this inequality might fail to be 
satisfied. First of all, if there is a localized region in the dielectric medium in 
which the refractive index suddenly changes {i.e., if there is an interface), then 
(4.180) is likely to break down in this region, allowing strong reflection of the 
incident wave. Secondly, the inequality obviously breaks down in the vicinity of a 
point where n = 0. We would, therefore, expect strong reflection of the incident 
wave from such a point. 



4.17 The reflection coefficient 

Consider an ionosphere in which the refractive index is a slowly varying function 
of height z above the surface of the Earth. Let be positive for z < zq, and 
negative for z > zq. Suppose that an upgoing radio wave of amplitude is 
generated at ground level {z = 0). The complex amplitude of the wave in the 
region < z < zq is given by the upgoing W.K.B. solution 

Ey = Eon-^/^exp ^ij^y ndz^ , (4.182a) 

cB^ = -Eon^/^exp ^ifc^ nd^j . (4.182b) 

The upgoing energy flux is given by -{EyB^* +Ey* B^) / (A^q) = (eo//io)^/^ l^^olV^- 
In the region z > zq the W.K.B. solutions take the form 



Ey = ^e*'^/^ |n|-^/^exp |^±fcy \n\dzj, (4.183a) 
cB^ = ±^e-''^/^|n|^/2exp ( ±jfe / \n\dz] , (4.183b) 
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where 71 is a constant. These solutions correspond to exponentially growing and 
decaying waves. Note that the magnetic components of the waves are in phase 
quadrature with the electric components. This implies that the Poynting fluxes 
of the waves are zero; i.e., the waves do not transmit energy. Thus, there is a 
non-zero incident energy flux in the region z < zq, and zero energy flux in the 
region z > zq. Clearly, the incident wave is either absorbed or reflected in the 
vicinity of the plane z = zq (where n = 0). In fact, as we shall prove later on, 
the wave is reflected. The complex amplitude of the reflected wave in the region 
< ^ < ^0 is given by the downgoing W.K.B. solution 

Ey = Eoi?n-^/^exp ^-i^y^ ndz^ , (4.184a) 

cB^ = EoRn^/'^exp(^-ik J^ndz^ , (4.184b) 

where R is the coeflicient of reflection. Suppose, for the sake of argument, that 
the plane z = zq acts like a perfect conductor, so that Ey{zo) = 0. It follows 
that 



\kjjndzy (4.185) 



R = — exp ^2i 
In fact, as we shall prove later on, the correct answer is 

R=-\ exp ( 2i / n dz] . (4.186) 



Thus, there is only a — 7r/2 phase shift at the reflection point, instead of the 
— TT phase shift which would be obtained if the plane z = zq acted like a perfect 
conductor. 



4.18 Extension to oblique incidence 

We have discussed the W.K.B. solutions for radio waves propagating vertically 
through an ionosphere whose refractive index varies slowly. Let us now generalize 
these solutions to allow for radio waves which propagate at an angle to the vertical 
axis. 
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The refractive index of the ionosphere varies continuously with height z. How- 
ever, let us, for the sake of clarity, imagine that the ionosphere is replaced by a 
number of thin discrete strata in which the medium is homogeneous. By mak- 
ing these strata sufficiently thin and numerous we can approximate as closely as 
is desired to the real ionosphere. Suppose that a plane wave is incident on the 
ionosphere, from below, and suppose that the wave normal lies in the x-z plane 
and makes an angle Qj with the vertical axis. At the lower boundary of the first 
stratum the wave is partially reflected and partially transmitted. The transmit- 
ted wave is partially reflected and partially transmitted at the second boundary 
between the strata, and so on. However, in the limit of many strata, where 
the diflFerence in refractive indices between neighbouring strata is very small, the 
amount of reflection at the boundaries becomes negligible. In the nth stratum, 
let n„ be the refractive index, and let Bn be the angle between the wave normal 
and the vertical axis. According to Snell's law, 

n„_i sin6'„_i = n„ sin6'„. (4.187) 

Below the ionosphere n = 1, and so 

Un sin^„ = sin ^7. (4.188) 

For a wave in the nth stratum, any field quantity depends on z and x through 
factors 

74exp [ i k n„(±2; cos ^„ + x sin ^„)] , (4.189) 

where ^4 is a constant. The ± signs denote upgoing and downgoing waves, re- 
spectively. When the operator djdx acts on the above expression, it is equivalent 
to multiplication by i/cn„ sin^„ = \ k sin^/, which is independent of x and z. It 
is convenient to use the notation S = sin^/. Hence, we may write symbolically 

d 

— = ikS, (4.190a) 
ox 

d 

— = 0. (4.190b) 
dy 

This result is true no matter how thin the strata are, so it must also hold for the 
real ionosphere. Note that, according to Snell's law, if the wave normal starts off 
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in the x-z plane then it will remain in this plane as it propagates through the 
ionosphere. 



Equations (4.164) and (4.190) can be combined to give 



and 



dcB^ 
dz 



dz 



= ikcB^, (4.191a) 
oz 

ikSEy = ikcB^, (4.191b) 

ikScB^ = -\kn^ Ey, (4.191c) 



-\kSEz = ikcBy, (4.192a) 

= -ikn'E,, (4.192b) 
ikScBy = -\kn^ E^. (4.192c) 



As before, Maxwell's equations can be split into two independent groups, corre- 
sponding to two independent polarizations of radio waves propagating through 
the ionosphere. For the first set of equations, the electric field is always parallel 
to the y-axis. The corresponding waves are, therefore, said to be horizontally 
polarized. For the second set of equations, the electric field always lies in the 
x-z plane. The corresponding waves are, therefore, said to be vertically polarized 
[n.b., the term "vertically polarized" does not necessarily imply that the electric 
field is parallel to the vertical axis). Note that the equations governing horizon- 
tally polarized waves are not isomorphic to those governing vertically polarized 
waves, so both types of waves must be dealt with separately. 

For the case of horizontally polarized waves, Eqs. (4.191b) and (4.191c) yield 

dcB, 



dz 

where 



= -ikq^Ey, (4.193) 
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The above equation can be combined with Eq. (4.191a) to give 



(4.195) 



Equations (4.193) and (4.195) have exactly the same form as Eqs. (4.166b) and 
(4.168), except that is replaced by , so the results of Section 4.16 can be 
immediately employed to find the W.K.B. solutions, which take the form 



Ey = Aq ^/^exp^±ifc J qdz^ , 



^Aq^^^ exp I ±ik qdz] , 



(4.196a) 



(4.196b) 



where A is a constant. Of course, both expressions should also contain a multi- 
plicative factor e^(^'^^~'^*), but this is usually omitted for the sake of clarity. By 
analogy with Eq. (4.180), the W.K.B. solutions are valid as long as 



1 



3 f]_dq 

4 \q'^ dz 



1 d^q 



2^3 dz'^ 



< 1. 



(4.197) 



This inequality clearly fails in the vicinity of ^ = 0, no matter how slowly q 
varies with z. Hence, q = 0, or = , specifies the height at which reflection 
takes place. By analogy with Eq. (4.186), the reflection coefficient at ground level 
{z = 0) is given by 

R= -iexp(2ik qdz^ , (4.198) 

where zq is the height at which q = 0. 

For the case of vertical polarization, Eqs. (4.192a) and (4.192c) yield 

dE, .,q' 
—— = ik — cBy. 
oz 

This equation can be combined with Eq. (4.192b) to give 

d'^By 1 d(n'^) dBy 2 2 



(4.199) 



dz' 



■n? dz dz 



+ fc^r By = 0. 



(4.200) 
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Clearly, the differential equation which governs the propagation of vertically po- 
larized waves is considerably more complicated than the corresponding equation 
for horizontally polarized waves. 

The W.K.B. solution for vertically polarized waves is obtained by substituting 
the wave-like solution 

cBy = Ae''^'^'\ (4.201) 

where ^ is a constant and (f){z) is the generalized phase, into Eq. (4.200). The 
differential equation thus obtained for the phase is 



.^20 /d(f)Y i d{n'^)d4> 



dz^ \dz) ■n? dz dz ^ ^ ^ ^ ^ ^ ^ 

Since the medium is slowly varying, the first and third term in the above equation 
are small, and so to a first approximation 

^ = ±kq, (4.203a) 
dz 

These expressions can be inserted into the first and third terms of Eq. (4.202) to 
give the second approximation 



dz 



\dz n dz J 



(4.204) 



The final two terms on the right-hand side of the above equation are small, so 
expanding the right-hand side using the binomial theorem yields 

d0 ±fc _^ ^ dq '\ dn ,^ 2Qb) 
dz 2q dz n dz 

This expression can be integrated, and the result inserted into Eq. (4.201), to 
give the W.K.B. solution 

cBy = Anq-^/"^ exp (^\k j qdz^ . (4.206) 
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The corresponding W.K.B. solution for is obtained from Eq. (4.199): 

= ±An-^q^l'^ exp ^±ifc ^ qdz^ . (4.207) 

Here, any terms involving derivatives of n and q have been neglected. 

Substituting Eq. (4.206) into the differential equation (4.200), and demanding 
that the result be small compared to the original terms in the differential equation, 
yields the following condition for the validity of the above W.K.B. solutions: 



1 




1 dq 



2q^ dz 



d^q 1 

+ ^ 



1 d'^n 



n dz' 



\n dz J 



< 1. 



(4.208) 



This criterion fails close to g' = 0, no matter how slowly n and q vary with z. 
Hence, g = gives the height at which reflection takes place. The condition also 
fails close to n = 0, which does not correspond to the reflection level. If, as is 
usually the case, the electron density in the ionosphere increases monotonically 
with height, then the level where n = lies above the reflection level, where q = {). 
If the two levels are well separated then the reflection process is unaffected by the 
failure of the above inequality at the level n = 0, and the reflection coefficient is 
given by Eq. (4.198), just as for the case of horizontal polarization. If, however, 
the level n = lies close to the level g = then the reflection coefficient may 
be affected, and a more accurate treatment of the differential equation (4.200) is 
required in order to obtain the true value of the reflection coefficient. 



4.19 Pulse propagation in the ionosphere 

Suppose that we possess a generator of radio waves which sends radio pulses 
vertically upwards into the ionosphere. For the sake of argument, we shall assume 
that these pulses are linearly polarized such that the electric field vector lies 
parallel to the y-axis. The pulse structure can be represented as 

/oo 
F(uj)q-'''^ duj, (4.209) 
-oo 
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where Ey{t) is the electric field produced by the generator [i.e., the field at z = 0). 
Suppose that the pulse is a signal of roughly constant (angular) frequency ujq, 
which lasts a time T, where T is long compared to I/ujq. It follows that F{uj) 
possesses narrow maxima around uj = ztujQ. In other words, only those frequencies 
which lie very close to the central frequency ujq play a significant role in the 
propagation of the pulse. 

Each component frequency of the pulse yields a wave which travels indepen- 
dently up into the ionosphere, in a manner specified by the appropriate W.K.B. 
solution [see Eqs. (4.181)]. Thus, if Eq. (4.209) specifies the signal at ground level 
(2; = 0), then the signal at height z is given by 

where 

Muj , z ,t) = — / n{uj, z) dz — ut. 
c Jo 

Here, we have used k = uj/c. 

Equation (4.210) can be regarded as a contour integral in a;-space. The quan- 
tity F/n^/^ is a relatively slowly varying function of cj, whereas the phase is a 
large and rapidly varying function of cj. As described in Section 4.11, the rapid 
oscillations of exp(i0) over most of the path of integration ensure that the in- 
tegrand averages almost to zero. However, this cancellation argument does not 
apply to those points on the path of integration where the phase is stationary] i.e., 
those points where d(j)/duj = 0. It follows that the left-hand side of Eq. (4.210) 
averages to a very small value, expect for those special values of z and t at which 
one of the points of stationary phase in w-space coincides with one of the peaks 
of F[uS). The locus of these special values of z and t can obviously be regarded 
as the equation of motion of the pulse as it propagates through the ionosphere. 
Thus, the equation of motion is specified by 



(4.211) 



(4.212) 
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which yields 



c Jo 



d{ujn) 
duj 



dz. 



(4.213) 



Suppose that the ^-velocity of a pulse of central frequency ujq at height z 
is given by Uz{ujo,z). The differential equation of motion of the pulse is then 
dt = dz/uz- This can be integrated, using the boundary condition z = at t = 0, 
to give the full equation of motion: 



t- r 

Jo Uz' 



(4.214) 



A comparison of Eqs. (4.213) and (4.214) yields 

d[uj n{uj 



Uz{ujQ,z) = c 



duj 



(4.215) 



a;=a;o 



The velocity Uz is usually called the group velocity. It is easily demonstrated that 
the above expression for the group velocity is entirely consistent with that given 
previously [see Eq. (4.135)]. 



The dispersion relation (4.164) yields 

n[u),z) = I 1 — 



1/2 



(4.216) 



in the limit where electron collisions are negligible. The phase velocity of radio 
waves of frequency uj propagating vertically through the ionosphere is given by 



Vz{u),z) = 



(uj,z) 



■1/2 



n 



(4.217) 



According to Eqs. (4.215) and (4.216), the corresponding group velocity is 



uz{uj,z) = c I l-^;j2- I 



1/2 



(4.218) 
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It follows that 



VzUz = C^- 



(4.219) 



Note that as the reflection point z = zq [defined as the solution of a; = ujp{zo)] is 
approached from below, the phase velocity tends to infinity, whereas the group 
velocity tends to zero. 

Let r be the time taken for the pulse to travel from the ground to the re- 
flection level, and back to the ground again. The product cr/2 is termed the 
equivalent height of reflection, and is denoted h{uj), since it is a function of the 
pulse frequency, u. The equivalent height is the height to which the pulse would 
have to go if it always traveled with the velocity c. Since we know that a pulse of 
dominant frequency u propagates at height z with the ^-velocity Uz{co, z) (this is 
true for both upgoing and downgoing pulses) , and also that the pulse is reflected 
at the height zq{lj), where cu = uJp{zo), it follows that 



Note that the equivalent height of reflection, h{uj), is always greater than the 
actual height of reflection, zo{uj), since the group velocity Uz is always less than 
the velocity of light. The above equation can be combined with Eq. (4.218) to 
give 



Note that the integrand diverges as the reflection point is approached, but the 
integral remains flnite. 

4.20 Determining the ionospheric electron density profile 

We can measure the equivalent height of the ionosphere in a fairly straightfor- 
ward manner, by timing how long it takes a radio pulse fired vertically upwards 




(4.220) 



Hence, 




(4.221) 




(4.222) 
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to return to ground level again. We can, therefore, determine the function h{uj) 
experimentally by performing this procedure many times over, using pulses of dif- 
ferent central frequencies. But, is it possible to use this information to determine 
the number density of free electrons in the ionosphere as a function of height? In 
mathematical terms, the problem is as follows. Does a knowledge of the function 



pZo (w) 

hiu) = / 
Jo 



CO 



(4.223) 



where Up^{zQ) = cj^, necessarily imply a knowledge of the function uj^{z)? Note, 
of course, that ujp^{z) oc N{z). 

Let up' = V and ujp{z) = u{z). Equation (4.223) then becomes 



^;-i/2/i(^i/2)= r 

JO 



[i;-w(z)]i/2' 



(4.224) 



where ^(2:0) = and u[z) < v for < z < zq- Let us multiply both sides of the 
above equation by {w — v)~^^'^/7r and integrate from v = to w. We obtain 



- v-'/^w - v)-'/^ h{v'/^) dv = - / 

Jo Jo Jo 



dz 



{w — t')i/2 [v — w)^/2 

(4.225) 

Consider the double integral on the right-hand side. The region of v-z space over 
which this integral is performed is sketched in Fig. 15. It can be seen that, as long 
as zq{v^^'^) is a monotonically increasing function of z, we can swap the order of 
integration to give 



TT 



dv 



'u{z) (w-vyi'^ {v-uyi'^ 



dz. 



(4.226) 



Here, we have used the fact that the curve z = zoiy^^'^) is identical with the curve 
V = u{z). Note that if Zq{v^/'^) is not a monotonically increasing function of v 
then we can still swap the order of integration, but the limits of integration are, 
in general, far more complicated than those indicated above. The integral over v 
in the above expression can be evaluated using standard methods (by making the 
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substitution v = w cos^ 6-\-u sin^ 6): the result is simply tt. Thus, the expression 
(4.226) reduces to zo{w'^/'^). It follows from Eq. (4.225) that 

zo{w'/') = - / V-'/' {w - v)-'/' h{v'/') dv. (4.227) 

Making the substitutions v = w sin^ a and w^^"^ = cu, we obtain 

2 /""/^ 

^o(a;) = - / h{uj sin a) da. (4.228) 

By definition, uj = ujp at the reflection level z = zo- Hence, the above equation 
reduces to 

2 f^/^ 

z{ujp) = — h{ijjp sin a) da. (4.229) 

7^ Jo 

Thus, we can obtain ^ as a function of ujp (and, hence, cUp as a function of 
z) simply by taking the appropriate integral of the experimentally determined 
function h{uj). Since ujp{z) oc ■\yN{z), this means that we can determine the 
electron number density proflle in the ionosphere provided we know the variation 
of the equivalent height of the ionosphere with pulse frequency. The constraint 
that zo{uj) must be a monotonically increasing function of uj translates to the 
constraint that N{z) must be a monotonically increasing function of z. Note that 
we can still determine N{z) from h{uj) for the case where the former function is 
non-monotonic, it is just a far more complicated procedure than that outlined 
above. Incidentally, the technique by which we have inverted Eq. (4.222), which 
specifies h{uj) as some integral over ujp{z), to give ujp{z) as some integral over 
h{uj) is known as Abel inversion. 



4.21 Ray tracing in the ionosphere 

Suppose that we possess a radio antenna which is capable of launching radio 
waves of constant frequency a; into the ionosphere at an angle to the vertical. Let 
us consider the paths traced out by these waves in the x-z plane. For the sake of 
simplicity, we shall assume that the waves are horizontally polarized, so that the 
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Figure 15: A sketch of the region of v-z space over which the integral on the 
right-hand side of Eq. (4.223) is evaluated 

electric field vector always lies parallel to the y-axis. The signal emitted by the 
antenna (located at ^ = 0) can be represented as 



the sake of clarity. Suppose that the signal emitted by the antenna is mostly 
concentrated in a direction making an angle 6i with the vertical. It follows that 
F[S) possesses a narrow maximum around S = Sq, where Sq = sinOj. 

If Eq. (4.230) represents the signal at ground level, then the signal at height 
z in the ionosphere is easily obtained by making use of the W.K.B. solutions for 
horizontally polarized waves described in Section 4.18. We obtain 




(4.230) 



where k = uj/c. Here, the e 



time dependence of the signal is neglected for 




(4.231) 



where 




(4.232) 
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Equation (4.231) is basically a line integral in 5'-space. The quantity Fjq^l'^ is a 
relatively slowly varying function of whereas the phase is a large and rapidly 
varying function of S. As described in Section 4.11, the rapid oscillations of 
exp(i0) over most of the path of integration ensure that the integrand averages 
almost to zero. In fact, only those points on the path of integration where the 
phase is stationary (i.e., where d4>/dS = 0) make a significant contribution to 
the integral. It follows that the left-hand side of Eq. (4.231) averages to a very 
small value, except for those special values of x and z at which one of the points 
of stationary phase in ^-space coincides with the peak of F{S). The locus of 
these special values of x and z can clearly be regarded as the trajectory of the 
radio signal emitted by the antenna as it passes through the ionosphere. Thus, 
the signal trajectory is specified by 



0, (4.233) 



which yields 



We can think of this equation as tracing the path of a ray of radio frequency 
radiation, emitted by the antenna at an angle 9i to the vertical (where Sq = 
sin^/), as it propagates through the ionosphere. 

Now 

q^ = n^- S\ (4.235) 
so the ray tracing equation becomes 

^ = ^ / c2 - (4.236) 

where S is the sine of the initial (i.e., at the antenna) angle of incidence of the ray 
with respect to the vertical axis. Of course, Eq. (4.236) only holds for upgoing 
rays. For downgoing rays, a simple variant of the previous analysis using the 
downgoing W.K.B. solutions yields 

x = S , + S / , (4.237) 
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where n{zo) = S. Thus, the ray ascends into the ionosphere after being launched 
from the antenna, reaches a maximum height zq above the surface of the Earth, 
and then starts to descend. The ray eventually intersects the Earth's surface 
again a horizontal distance 

away from the antenna. 

The angle ^ which the ray makes with the vertical is given by tan^ = dx/dz. 
It follows from Eqs. (4.236) and (4.237) that 

S 

tan^ = ±— (4.239) 

where the upper and lower signs correspond to the upgoing and downgoing parts 
of the ray trajectory, respectively. Note that ^ = tt/2 at the reflection point, 
where n = S. Thus, the ray is horizontal at the reflection point. 

Let us investigate the reflection process in more detail. In particular, we wish 
to prove that radio waves are reflected at the ^ = surface, rather than being 
absorbed. We would also like to understand the origin of the — 7r/2 phase shift 
of radio waves at reflection which is evident in Eq. (4.198). In order to achieve 
these goals, we shall need to review the mathematics of asymptotic series. 



4.22 Asymptotic series: A mathematical aside 



It is often convenient to expand a function of the complex variable f{z) in inverse 
powers of z: 



^0 + — + ^ + 
z z^ 



(4.240) 



where 0(z) is a function whose behaviour for large values of z is known. If 
f{z)/(j){z) is singular as \z\ — > oo then the above series diverges. Nevertheless, 
under certain circumstances, the series may still be useful. 
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The circumstance needed to make this possible is that the difference between 
f{z)/(f){z) and the first n + 1 terms of the series be of order 1/^"+^, so that for 
sufficiently large z this difference becomes vanishingly small. More precisely, the 
series is said to represent f{z)/(f){z) asymptotically, that is 



p=0 



(4.241) 



provided that 



lim 



z 



n 



m 

<l>{z) 



A 

ZP 



p=0 



^0. 



(4.242) 



In other words, for a given value of n, the first n + 1 terms of the series may be 
made as close as desired to the ratio f{z)/4>{z) by making z sufficiently large. For 
each value of z and n there is an error of order 1 / z'^'^'^ . Since the series actually 
diverges, there is an optimum number of terms in the series used to represent 
f{z)/(f){z) for a given value of z. Associated with this is an unavoidable error. As 
z increases, the optimal number of terms increases and the error decreases. 



Consider a simple example. The exponential integral is defined 



El 



J X 



dt. 



(4.243) 



The asymptotic series for this function can be generated via a series of partial 
integrations. For example. 



El (a;) = 



X 



t2 



dt. 



(4.244) 



Continuing this procedure yields 



Ei(x) 



X 



1 2! 3! 



+ 



(-l)"n! 



X' 



r°° e-* 

+ (_!)"+!(„ +1),/ 



dt. 



(4.245) 
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The infinite series obtained by taking the hmit n — > oo diverges, since the Cauchy 
convergence test yields 



hm 


Un+1 


= hm 


n' 


— > oo 


n—^oo 


Un 


n—^oo 


■ X. 





(4.246) 



Note that two successive terms in the series become equal in magnitude for n = 
indicating that the optimum number of terms for a given x is roughly the integer 
nearest x. To prove that the series is asymptotic, we need to show that 



lim a;"+ie^(-l)"+i(n + l)! / 

Jx 

This immediately follows, since 

1 



oo g-t 



in+2 



dt = 0. 



(4.247) 



J^x 



oo g-t 



dt < 



X 



n+2 



dt 



X 



n+2 



(4.248) 



Thus, the error involved in using the first n terms is less than (n + 1)! e~^/a:""'"^, 
which is exactly the next term in the series. We can see that as n increases, this 
estimate of the error first decreases and then increases without limit. In order to 
visualize this phenomenon more exactly, let f{x) = a:;exp(a:;) ^{x), and let 



p=0 



(4.249) 



be the asymptotic series representation of this function which contains n+1 terms. 
Figure 16 shows the relative error in the asymptotic series \fn{x) — f{x)\/f{x) 
plotted as a function of the approximate number of terms in the series n for 
X = 10. It can be seen that as n increases the error initially falls, reaches a 
minimum value at about n = 10, and then increases rapidly. Clearly, the optimum 
number of terms in the asymptotic series used to represent /(lO) is about 10. 

Asymptotic series are fundamentally different to conventional power law ex- 
pansions, such as 



^3 ^3 ^7 

sin. = .-- + --- + 



(4.250) 
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Figure 16: The relative error in a typical asymptotic series plotted as a function 
of the number of terms in the series 



This series representation of sin^ converges absolutely for all finite values of z. 
Thus, at any z the error associated with the series can be made as small as is 
desired by including a sufficiently large number of terms. In other words, unlike 
an asymptotic series, there is no intrinsic, or unavoidable, error associated with 
a convergent series. It follows that a convergent power law series representation 
of a function is unique inside the domain of convergence of the series. On the 
other hand, an asymptotic series representation of a function is not unique. It is 
perfectly possible to have two different asymptotic series representations of the 
same function, as long as the difference between the two series is less than the 
intrinsic error associated with each series. Furthermore, it is often the case that 
different asymptotic series are used to represent the same single-valued analytic 
function in different regions of the complex plane. 

For example, consider the asymptotic expansion of the confiuent hypergeomet- 
ric function F{a,c,z). This function is the solution of the differential equation 

zF" + (c - z)F' -aF = (4.251) 

which is analytic at z = [in fact, F{a,c,0) = 1]. Here, ' denotes d/dz. The 
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asymptotic expansion of F{a, c, z) takes the form: 



r(a)r(c a)^^^^^^^^ ^ r(c-a)z"-'^eqi + 0(l/z)] 



r(c) 

for — TT < arg(^) < 0, and 
T{a)r(c - a) 



+r(a) e-''^" [1 + 0(l/z)] (4.252a) 



i^(a,c,z) ~ r(c-a)^«-^eqi + 0(l/z)] 

r(c) 



+r(a) e'"" [1 + 0(l/z)] (4.252b) 



for < arg(^) < tt, and 



r(a)r(c X ^ r(c-a)^"-^e-'2"("-^)eqi + 0(l/z)] 

r(c) 

+r(a) e'"" [1 + 0{l/z)] (4.252c) 

for TT < arg(2;) < 27r, etc. It can be seen that the expansion consists of a Hnear 
combination of two asymptotic series (only the first term in each series is shown). 
For \z\ 1, the first series is exponentially larger than the second whenever 
Ke{z) > 0. We say that the first series is dominant in this region, whereas the 
second series is subdominant. Likewise, the first series is exponentially smaller 
than the second whenever Re(z) < 0. We say that the first series is subdominant 
and the second series is dominant in this region. 

Consider a region in which one or other of the series is dominant. Strictly 
speaking, it is not mathematically consistent to include the subdominant series 
in the asymptotic expansion because its contribution is actually less than the 
intrinsic error associated with the dominant series [this error is 0{l/z) times the 
dominant series, since we are only including the first term in this series]. Thus, 
at a general point in the complex plane the asymptotic expansion simply consists 
of the dominant series. However, this is not the case in the immediate vicinity 
of the lines Re(^) = 0: these are called the anti- Stokes lines. When an anti- 
Stokes line is crossed, a dominant series becomes subdominant and vice versa. In 
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the immediate vicinity of an anti- Stokes line neither series is dominant, so it is 
mathematically consistent to include both series in the asymptotic expansion. 

The hypergeometric function F{a, c, z) is a perfectly well behaved, single- 
valued, analytic function in the complex plane. However, our two asymptotic se- 
ries are, in general, multi- valued functions in the complex plane [see Eq. (4.252a)]. 
Can a single-valued function be represented asymptotically by a multi-valued 
function? The short answer is no. We have to employ different combinations of 
the two series in different regions of the complex plane in order to ensure that 
F{a, c, z) remains single- valued. Equations (4.252) show how this is achieved. Ba- 
sically, the coefficient in front of the subdominant series changes dis continuously 
at certain values of arg(2;). This is perfectly consistent with F{a,c,z) being an 
analytic function because the subdominant series is "invisible"; i.e., the contribu- 
tion of the subdominant series to the asymptotic solution falls below the intrinsic 
error associated with the dominant series, so it does not really matter if the co- 
efficient in front of the former series changes discontinuously. Imagine tracing a 
large circle, centred on the origin, in the complex plane. Close to an anti- Stokes 
line, neither series is dominant, so we must include both series in the asymptotic 
expansion. As we move away from the anti-Stokes line, one series becomes dom- 
inant, which means that the other series becomes subdominant and, therefore, 
drops out of our asymptotic expansion. Eventually, we approach a second anti- 
Stokes line, and the subdominant series reappears in our asymptotic expansion. 
However, the coefficient in front of the subdominant series when it reappears is 
different to that which it had when it disappeared. This new coefficient is carried 
across the second anti-Stokes line into the region where the subdominant series 
becomes dominant. In this new region, the dominant series becomes subdominant 
and disappears from our asymptotic expansion. Eventually, a third anti-Stokes 
line is approached and the series reappears, but, again, with a different coefficient 
in front. The jumps in the coefficients of the subdominant series are chosen in 
such a manner that if we perform a complete circuit in the complex plane then 
the value of the asymptotic expansion is the same at the beginning and the end 
points. In other words, the asymptotic expansion is single- valued, despite the 
fact that it is built up out of two asymptotic series which are not single- valued. 
The jumps in the coefficient of the subdominant series, which are needed to keep 
the asymptotic expansion single-valued, are called Stokes phenomena, after the 
celebrated nineteenth century British mathematician Sir George Gabriel Stokes, 
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who first drew attention to this effect. 

Where exactly does the jump in the coefficient of the subdominant series 
occur? All we can really say is "somewhere in the region between two anti-Stokes 
lines where the series in question is subdominant." The problem is that we only 
retain the first term in each asymptotic series. Consequently, the intrinsic error 
in the dominant series is relatively large and we lose track of the subdominant 
series very quickly after moving away from an anti-Stokes line. However, we 
could include more terms in each asymptotic series. This would enable us to 
reduce the intrinsic error in the dominant series and, thereby, expand the region 
of the complex plane in the vicinity of the anti-Stokes lines where we can see 
both the dominant and subdominant series. If we were to keep adding terms 
to our asymptotic series, so as to minimize the error in the dominant solution, 
we would eventually be forced to conclude that a jump in the coefficient of the 
subdominant series can only take place on those lines in the complex plane on 
which Im(^) = 0: these are called Stokes lines. This result was first proved by 
Stokes in 1857.^^ On a Stokes line the magnitude of the dominant series achieves 
its maximum value with respect to that of the subdominant series. Once we know 
that a jump in the coefficient of the subdominant series can only take place at a 
Stokes line, we can retain the subdominant series in our asymptotic expansion in 
all regions of the complex plane. What we are basically saying is that, although, 
in practice, we cannot actually see the subdominant series very far away from an 
anti-Stokes line because we are only retaining the first term in each asymptotic 
series, we could, in principle, see the subdominant series at all values of arg(^) 
provided that we retained a sufficient number of terms in our asymptotic series. 

Figure 17 shows the location in the complex plane of the Stokes and anti- 
Stokes lines for the asymptotic expansion of the hypergeometric function. Also 
shown is a branch cut, which is needed to make z single- valued. The branch cut 
is chosen such that arg(^) = on the positive real axis. Every time we cross an 
anti-Stokes line the dominant series becomes subdominant and vice versa. Every 
time we cross a Stokes line the coefficient in front of the dominant series stays 
the same, but that in front of the subdominant series jumps discontinuously [see 
Eqs. (4.252)]. Finally, the jumps in the coefficient of the subdominant series are 
such as to ensure that the asymptotic expansion is single- valued. 

i^G.G. Stokes, Trans. Camb. Phil. Soc. 10, 106-128 (1857) 
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anti-Stokes line 




Figure 17: The location of the Stokes hnes (dashed), the anti-Stokes Hnes (soHd), 
and the branch cut (wavy) in the complex plane for the asymptotic expansion of 
the hypergeometric function 
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4.23 The W.K.B. solutions as asymptotic series 



We have seen that the W.K.B. solution 




(4.253) 



(4.254) 



in the limit where the typical wavelength, 2'K/nk, is much smaller than the typical 
variation length-scale of the refractive index. But, what sort of approximation is 
involved in writing this solution? 

It is convenient to define the scaled variable 



where L is the typical variation length-scale of the refractive index, n{z). Equa- 
tion (4.254) can then be written 



where w{z,h) = Ey{Lz), q{z) = n'^{Lz), ' = d/dz^ and h = kL. Note that, in 
general, q{z), q"{z), etc. are 0(1) quantities. The non-dimensional constant 

h is of order the ratio of the variation length-scale of the refractive index to the 
wavelength. Let us seek the solutions to Eq. (4.256) in the limit 1. 

We can write 




z 



(4.255) 



w" -\- qw = 0, 



(4.256) 



w{z, h) = 
Equation (4.256) transforms to 



exp [ i h (f){z, h)] . 



(4.257) 



Expanding in powers of l//i, we obtain 




(0r +^=0. 



(4.258) 




(4.259) 
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which yields 



w 



{z,h) 



■1/4 



exp 




(4.260) 



Of course, we immediately recognize this expression as a W.K.B. solution. 

Suppose that we keep expanding in powers of 1/h in Eq. (4.259). The appro- 
priate generalization of Eq. (4.260) is a series solution of the form 



w 



{z,h) 



-1/4 

q ' exp 



±i 



qdz 



oo 



1 + E 

p=i 



hP 



(4.261) 



This is, in fact, an asymptotic series in h. We can now appreciate that a W.K.B. 
solution is just a highly truncated asymptotic series in h, in which only the first 
term in the series is retained. 

But, why is it so important that we recognize that W.K.B. solutions are 
highly truncated asymptotic series? The point is that the W.K.B. method was 
initially rather controversial after it was popularized in the 1920s. A lot of people 
thought that the method was completely wrong. Let us try to understand what 
the problem was. Suppose that we have never heard of an asymptotic series. 
Looking at Eq. (4.261), we would imagine that the expression in square brackets 
is a power law expansion in 1/h. The W.K.B. approximation involves neglecting 
all terms in this expansion except the first. This sounds fine, as long as h is 
much greater than unity. But, surely, to be mathematically rigorous, we have to 
check that the sum of all of the terms in the expansion which we are neglecting 
is small compared to the first term? However, if we attempt this we discover, 
much to our consternation, that the expansion is divergent. In other words, the 
sum of all of the neglected terms is infinite! Thus, if we interpret Eq. (4.261) 
as a conventional power law expansion in 1/h, the W.K.B. method is clearly 
nonsense: the W.K.B. solution is the first approximation to infinity. However, 
once we appreciate that Eq. (4.261) is actually an asymptotic series in h, the 
fact that the series diverges becomes irrelevant. If we retain the first n terms 
in the series, the series approximates the exact solution of Eq. (4.261) with an 
intrinsic (fractional) error which is of order l//i" {i.e., the first neglected term in 
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the series). The error is minimized at a particular value of h. As the number of 
terms in the series is increased, the intrinsic error decreases, and the value of h 
at which the error is minimized increases. In particular, we can see that there is 
an intrinsic error associated with a W.K.B. solution which is of order 1/h times 
the solution. 

It is amusing to note that if Eq. (4.261) were not a divergent series then it 
would be impossible to obtain total reflection of the W.K.B. solutions at the point 
^ = 0. As we shall discover, the reflection is directly associated with the fact that 
the expansion (4.261) exhibits a Stokes phenomenon. It is, of course, impossible 
for a convergent power series expansion to exhibit a Stokes phenomenon. 



4.24 Stokes constants 



We have seen that the differential equation 

w" + h^q(z)w = 0, 
where ' = d/dz, possesses approximate W.K.B. solutions of the form 



(4.262) 



{a,z) = q exp \ ih q^^^ dz 



1 + 



h 



(i,a) 



1 /., 

q ' exp 



Ah 



^1/2 dz 



1 + 



h 



(4.263a) 



(4.263b) 



Here, we have adopted an arbitrary phase reference level z = a. The convenient 
notation (a, z) is fairly self explanatory: a and z refer to the lower and upper 
bounds of integration, respectively, inside the exponential. It follows that the 
other W.K.B. solution can be written (i,a) (we can reverse the limits of inte- 
gration inside the exponential to obtain minus the integral in z from z = a to 
z = z). 

Up to now we have thought of z as a real variable representing scaled height 
in the ionosphere. Let us now generalize our analysis somewhat and think of z 
as a CO mp /ex variable. There is nothing in our derivation of the W.K.B. solutions 
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which depends crucially on z being a real variable, so we expect these solutions to 
remain valid when z is reinterpreted as a complex variable. Incidentally, we must 
now interpret q{z) as some well behaved function of the complex variable. An 
approximate general solution of the differential equation (4.262) in the complex 
z plane can be written linear superposition of the two W.K.B. solutions 

(4.263). 

The parameter h is assumed to be much larger than unity. It is clear from 
Eqs. (4.263) that in some regions of the complex plane one of the W.K.B. solutions 
is going to be exponentially larger than the other. In such regions, it is not 
mathematically consistent to retain the smaller W.K.B. solution in the expression 
for the general solution, since the contribution of the smaller W.K.B. solution is 
less than the intrinsic error associated with the larger solution. Adopting the 
terminology introduced in Section 4.22, the larger W.K.B. solution is said to be 
dominant, and the smaller solution is said to be subdominant Let us denote 
the W.K.B. solution (4.263a) as {a,z)d in regions of the complex plane where it 
is dominant, and as (a, i)s in regions where it is subdominant. An analogous 
notation is adopted for the second W.K.B. solution (4.263b). 

Suppose that q{z) possesses a simple zero at the point z = zq (chosen to 
be the origin for the sake of convenience). It follows that in the immediate 
neighbourhood of the origin we can write 

q = aiz-\-a2z'^ -\----, (4.264) 

where ai ^ 0. It is convenient to adopt the origin as the phase reference point 
{i.e., a = 0), so the two W.K.B. solutions become (0, z) and (i, 0). We can define 
anti-Stokes lines in the complex z plane (see Section 4.22). These are lines which 
satisfy 

= 0. (4.265) 

As we cross an anti-Stokes line, a dominant W.K.B. solution becomes subdom- 
inant, and vice versa. Thus, (0,i)c^ ^ (0,i)s and {z,0)ci ^ (i,0)s. In the 
immediate vicinity of an anti-Stokes line the two W.K.B. solutions have about 
the same magnitude, so it is mathematically consistent to include the contribu- 
tions from both solutions in the expression for the general solution. In such a 
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region, we can drop the subscripts d and s, since the W.K.B. solutions are neither 
dominant nor subdominant, and write the W.K.B. solutions simply as (0, z) and 
(^,0). 

It is clear from Eqs. (4.263) that the W.K.B. solutions are not single-valued 
functions of i, since they depend on g^/^(i), which is a double- valued function. 
Thus, if we wish to write an approximate analytic solution to the differential 
equation (4.262) we cannot express this solution as the same linear combination 
of W.K.B. solutions in all regions of the complex ^-plane. This implies that there 
must exist certain lines in the complex i-plane across which the mix of W.K.B. 
solutions in our expression for the general solution changes discontinuously. These 
lines are called Stokes lines (see Section 4.22), and satisfy 

= 0. (4.266) 

As we cross a Stokes line, the coefficient of the dominant W.K.B. solution in our 
expression for the general solution must remain unchanged, but the coefficient of 
the subdominant solution is allowed to change discontinuously. Incidentally, this 
is perfectly consistent with the fact that the general solution is analytic: the jump 
in our expression for the general solution due to the jump in the coefficient of the 
subdominant W.K.B. solution is less than the intrinsic error in this expression 
due to the intrinsic error in the dominant W.K.B. solution. Once we appreciate 
that the coefficient of the subdominant solution can only change at a Stokes line, 
we can retain both W.K.B. solutions in our expression for the general solution 
throughout the complex z plane. In practice, we can only see a subdominant 
solution in the immediate vicinity of an anti- Stokes line, but if we were to evaluate 
the W.K.B. solutions to higher accuracy [i.e. retain more terms in the asymptotic 
series in Eqs. (4.263)]we could, in principle, follow a subdominant solution all the 
way to a neighbouring Stokes line. 

In the immediate vicinity of the origin 

qU^dzc^'^z'/'. (4.267) 

It follows from Eqs. (4.265) and (4.266) that three Stokes lines and three anti- 
Stokes lines radiate from a zero of q{z). The general arrangement of Stokes and 
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Figure 18: The arrangement of Stokes lines (dashed) and anti-Stokes hnes (sohd) 
around a simple zero of q{z). Also shown is the branch cut (wavy line). All of 
the lines radiate from the point q = 0. 



anti-Stokes lines in the vicinity of a g = point is sketched in Fig. 18. Note that 
a branch cut must also radiate from the q = point in order to uniquely specify 
the function g^/^(i). Thus, in general, seven lines radiate from a zero of q{z), 
dividing the complex z plane into seven domains (numbered 1 through 7) . 

Let us write our general solution as 

w{z, h) = A (0, z)+B (i, 0) (4.268) 

on the anti-Stokes line between domains 1 and 7, where A and B are arbitrary 
constants. Suppose that the W.K.B. solution (0,i) is dominant in domain 7. 
Thus, in domain 7 the general solution takes the form 

w{7)=A{0,z)d + B{z,0)s. (4.269) 

Let us move into domain 1. In doing so, we cross an anti-Stokes line, so the 
dominant solution becomes subdominant, and vice versa. Thus, in domain 1 the 
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general solution takes the form 



w{l) = A (0, z)s + B {z, 0)d. (4.270) 

Let us now move into domain 2. In doing so, we cross a Stokes line, so the 
coefficient of the dominant solution, 5, must remain constant, but the coefficient 
of the subdominant solution. A, is allowed to change. Suppose that the coefficient 
of the subdominant solution jumps by t times the coefficient of the dominant 
solution, where t is an undetermined constant. It follows that in domain 2 the 
general solution takes the form 

w(2) = (A + t B) (0, z), + B (i, 0)d. (4.271) 

Let us now move into domain 3. In doing so, we cross an anti-Stokes line, so the 
the dominant solution becomes subdominant, and vice versa. Thus, in domain 3 
the general solution takes the form 

w;(3) = (A + tB)(0,i)d + B(i,0),. (4.272) 

Let us now move into domain 4. In doing so, we cross a Stokes line, so the 
coefficient of the dominant solution must remain constant, but the coefficient of 
the subdominant solution is allowed to change. Suppose that the coefficient of the 
subdominant solution jumps by u times the coefficient of the dominant solution, 
where u is an undetermined constant. It follows that in domain 4 the general 
solution takes the form 

w{A) = {A + tB){0,z)d + {B^u[A + tB]){z,0)s. (4.273) 

Let us now move into domain 5. In doing so, we cross an anti-Stokes line, so the 
the dominant solution becomes subdominant, and vice versa. Thus, in domain 5 
the general solution takes the form 

w{5) = (A + t B) {{),z)s + (B + u[A + tB]) (i, 0)^. (4.274) 

Let us now move into domain 6. In doing so, we cross the branch cut in an anti- 
clockwise direction. Thus, the argument of z decreases by 27r. It follows from 
Eq. (4.264) that ^^/^ —q^l'^ and q^l^ —'\q^l^. The following rules for tracing 
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the W.K.B. solutions across the branch cut (in an anti-clockwise direction) ensure 
that the general solution is continuous across the cut [see Eqs. (4.261)]: 



(o,i) ^ -i(i,o), 
(i,o) ^ -i(o,i). 



(4.275a) 
(4.275b) 



Note that the properties of dominancy and subdominancy are preserved when 
the branch cut is crossed. It follows that in domain 6 the general solution takes 
the form 



Let us, finally, move into domain 7. In doing so, we cross a Stokes line, so the 
coefficient of the dominant solution must remain constant, but the coefficient of 
the subdominant solution is allowed to change. Suppose that the coefficient of the 
subdominant solution jumps by v times the coefficient of the dominant solution, 
where v is an undetermined constant. It follows that in domain 7 the general 
solution takes the form 



w{l) = -i{A+tB + v{B + u [A+tB]}) (i, 0), -i {B+u [A+t B]) (0, z)^. (4.277) 



Now, we expect our general solution to be an analytic function, so it follows 
that the solutions (4.269) and (4.277) must be identical. Thus, we can compare 
the coefficients of the two W.K.B. solutions, (i, 0)s and (0,i)d. Since A and B 
are arbitrary, we can also compare the coefficients of A and B. Thus, comparing 
the coefficients of A (0, z)d, we find 



w;(6) 



i{A^tB) (i, Q),-i{B + u[A + tB]) (0, z)d. 



(4.276) 



1 = -i 



(4.278) 



Comparing the coefficients of B (0, z)^ yields 



= 1 + ut. 



(4.279) 



Comparing the coefficients of 71 (i, 0)s gives 



(4.280) 
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Finally, comparing the coefficients of B {z,0)s yields 
Equations (4.278) (4.281) imply that 



(4.281) 



t = u = V = i. 



(4.282) 



In other words, if we adopt the simple rule that every time we cross a Stokes line 
in an anti-clockwise direction the coefficient of the subdominant solution jumps 
by i times the coefficient of the dominant solution, then this ensures that our 
expression for the general solution (4.268) behaves as an analytic function. Here, 
the constant i is usually called a Stokes constant. Note that if we cross a Stokes 
line in a clockwise direction then the coefficient of the subdominant solution has 
to jump by — i times the coefficient of the dominant solution in order to ensure 
that our general solution behaves as an analytic function. 

4.25 The reflection coefficient 

Let us write z = x -\- iy, where x and y are real variables. Consider the solution 
of the differential equation 



where q{x) is a real function, /i is a large number, q > for x < 0, and ^ < for 
X > 0. It is clear that z = represents a simple zero of q{z). Here, we assume, 
as seems eminently reasonable, that we can find a well behaved function of the 
complex variable q{z) such that q{z) = q{x) along the real axis. The arrangement 
of Stokes and anti- Stokes lines in the immediate vicinity of the point i = is 
sketched in Fig. 19. The argument of q{z) on the positive x-axis is chosen to be 
— TT. Thus, the argument of q{z) on the negative x-axis is 0. 

On OB, the two W.K.B. solutions (4.261) can be written 



w" + q{x) w = 0, 



(4.283) 




(4.284a) 




(4.284b) 
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Figure 19: The arrangement of Stokes lines (dashed) and anti-Stokes hnes (sohd) 
in the complex z plane. Also shown is the branch cut (wavy line). 



206 



Here, we can interpret (0, x) as a wave propagating to the right along the x-axis, 
and {x, 0) as a wave propagating to the left. On OA, the W.K.B. solutions take 
the form 




(4.285a) 
(4.285b) 



Clearly, {x,0)s represents an evanescent wave which decays to the right along the 
a;-axis, whereas (O,^)^^ represents an evanescent wave which decays to the left. If 
we adopt the boundary condition that there is no incident wave from the region 
X — >• +00, the most general asymptotic solution to Eq. (4.283) on OA is written 

w{x,h) =A(x,0)s, (4.286) 

where A is an arbitrary constant. 

Let us assume that we can find an analytic solution w{z,h) to the diflFerential 
equation 

w" + q{z)w = {}, (4.287) 

which satisfies w{z, h) = w{x, h) along the real axis, where w{x, h) is the physical 
solution. From a mathematical point of view, this seems eminently reasonable. 
In the domains 1 and 2 the solution (4.286) becomes 

w{l) = A{z,0)s, (4.288) 

and 

w{2) = A{z,0)s. (4.289) 

Note that the solution is continuous across the Stokes line OA, since the coeflficient 
of the dominant solution (0,i) is zero: thus, the jump in the coefficient of the 
subdominant solution is zero times the Stokes constant, i; i.e., it is zero. Let us 
move into domain 3. In doing so, we cross an anti-Stokes line, so the solution 
becomes 

w{3) = A{z,0)d. (4.290) 
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Let us now move into domain 4. In doing so, we cross a Stokes line. Applying 
the general rule derived in the preceding section, the solution becomes 

w{4) = A{z,0)d^iA{0,z),. (4.291) 

Finally, on OB the solution becomes 

w{x,h) = ^(a;,0) + i^(0,a;). (4.292) 



Suppose that there is a point a on the negative a:- axis where q{x) = 1. It 
follows from Eqs. (4.286) and (4.292) that we can write the asymptotic solution 
to Eq. (4.283) as 

w{x,h) = g-^/^ exp l^i/^y^ q^/'^dx^ (4.293) 

-iexp^2i/i^ q^/'^dx^ ^-^/^ exp ^-i /i q^/'^ dx^ , 
in the region x < 0, and 



w 



{x,h)=e^p(^ih q^/'^dx^ e"''^/^ exp ^-/i ^ \q\^^'^ dx^ (4.294) 

in the region x > 0. Here, we have chosen 

^ = -i exp ^ i /i J ^1/2 dx^ . (4.295) 

If we interpret normalized altitude in the ionosphere, q{x) as the square 

of the refractive index in the ionosphere, the point a as ground level, and w as 
the electric field strength of a radio wave propagating vertically upwards into the 
ionosphere, then Eq. (4.293) tells us that a unit amplitude wave fired vertically 
upwards from ground level into the ionosphere is reflected at the level where the 
refractive index is zero. The first term in Eq. (4.293) is the incident wave and 
the second term is the reflected wave. The reflection coefl&cient {i.e., the ratio of 
the reflected to the incident wave at ground level) is given by 

= -i exp ^2 i /i J ^1/2 dx^ . (4.296) 
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Note that \R\ = 1, so the ampHtude of the reflected wave equals that of the 
incident wave. In other words, there is no absorption of the wave at the level 
of reflection. The phase shift of the reflected wave at ground level, with respect 
to that of the incident wave, is that associated with the wave propagating from 
ground level to the reflection level and back to ground level again, plus a — 7r/2 
phase shift at reflection. According to Eq. (4.294), the wave attenuates fairly 
rapidly (in the space of a few wavelengths) above the reflection level. Of course, 
Eq. (4.296) is completely equivalent to Eq. (4.186). 

Note that the reflection of the incident wave at the point where the refractive 
index is zero is directly associated with the Stokes phenomenon. Without the 
jump in the coefiicient of the subdominant solution, as we go from domain 3 to 
domain 4, there is no reflected wave on the OB axis. Note, also, that the W.K.B. 
solutions (4.293) and (4.294) break down in the immediate vicinity of q = [i.e., 
the reflection point). Thus, it is possible to demonstrate that the incident wave 
is totally reflected at the point ^ = 0, with a — 7r/2 phase shift, without having 
to solve for the wave structure in the immediate vicinity of the reflection point. 
This demonstrates that the reflection of the incident wave at ^ = is an intrinsic 
property of the W.K.B. solutions, and does not depend on the detailed behaviour 
of the wave in the region where the W.K.B. solutions break down. 

4.26 The Jeffries connection formula 

In the preceding section there is a tacit assumption that the square of the refrac- 
tive index, q{x) = n^(x), is a real function. As is apparent from Eq. (4.162), this 
is only the case in the ionosphere as long as electron collisions are negligible. Let 
us generalize our analysis to take electron collisions into account. In fact, the 
main effect of electron collisions is to move the zero of q{z) a short distance off 
the real axis (the distance is relatively short provided that we adopt the physical 
ordering z/ <C a;). The arrangement of Stokes and anti-Stokes lines around the 
new zero point, located at ^ = is sketched in Fig. 20. Note that electron 
collisions only significantly modify the form of q{z) in the immediate vicinity of 
the zero point. Thus, sufficiently far away from z = zq in the complex i-plane, 
the W.K.B. solutions and the locations of the Stokes and anti-Stokes lines are 
exactly the same as in the preceding section. 
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Figure 20: The arrangement of Stokes lines (dashed) and anti- Stokes hnes (sohd) 
in the complex z plane. Also shown is the branch cut (wavy line). 
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The W.K.B. solutions (4.284) and (4.285) are valid all the way along the 
real axis, except for a small region close to the origin where electron collisions 
significantly modify the form of q{z). Thus, we can still adopt the physically 
reasonable decaying solution (4.286) on the positive real axis. Let us trace this 
solution in the complex i-plane until we reach the negative real axis. We can 
achieve this by moving in a semi- circle in the upper half-plane. Since we never 
move out of the region in which the W.K.B. solutions (4.284) and (4.285) are 
valid, we conclude, by analogy with the preceding section, that the solution on 
the negative real axis is given by Eq. (4.292). Of course, in all of the W.K.B. 
solutions the point i = must be replaced by the new zero point z = zq. The new 
formula for the reflection coefficient, which is just a straightforward generalization 



This is called the Jeffries connection formula, after H. Jeffries, who discovered it 
in 1923. The general expression for the reflection coefficient is incredibly simple. 
We just integrate the W.K.B. solution in the complex i-plane from the phase 
reference level z = a to the zero point, square the result, and multiply by — i. 
Note that the path of integration between z = a and z = zq does not matter, 
because of Cauchy's theorem. Note, also, that since q^^^ is, in general, complex 
along the path of integration, we no longer have \R\ = 1. In fact, it is easily 
demonstrated that \R\ < 1. Thus, when electron collisions are included in the 
analysis we no longer obtain perfect reflection of radio waves from the ionosphere. 
Instead, some (small) fraction of the radio energy is absorbed at each reflection 
event. This energy is ultimately transfered to the particles in the ionosphere with 
which the electrons collide. 



of Eq. (4.296), is 




(4.297) 
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